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TRANSPORT  PHENOMENA  AND  INTERFACIAL  KINETICS 
IN  MULTIPHASE  COMBUSTION  SYSTEMS 

Principal  Investigator:  Prof.  Daniel  E.  Rosner 
1.  INTRODUCTION 

The  performance  of  ramjets  burning  slurry  fuels  (leading  to  condensed  oxide  aerosols  and 
liquid  film  deposits),  gas  turbine  engines  in  dusty  or  marine  atmospheres,  or  when  using  fuels 
from  non-traditional  sources,  depends  upon  the  formation  and  transport  of  small  particles  across 
non-isothermal  combustion  gas  boundary  layers  (BLs).  Even  airbreathing  engines  burning 
"clean"  hydrocarbon  fuels  can  experience  5001  formation/deposition  problems  (e.g.,  combustor 
liner  burnout,  accelerated  turbine  blade  erosion  and  "hot"  corrosion).  Moreover,  particle 
formation  and  transport  are  important  in  many  chemical  reactors  used  to  synthesize  or  process 
aerospace  materials  (turbine  blade  coatings,  optical  waveguides,  ceramic  precursor  powders, 
fibers  for  composites,...).  Accordingly,  our  research  is  directed  toward  providing  chemical 
propulsion  systems  engineers  and  materials-oriented  engineers  with  new  techniques  and 
quantitative  information  on  important  particle-  and  vapor-mass  transport  mechanisms  and  rates. 

The  purpose  of  this  report  is  to  summarize  our  research  methods  and  accomplishments 
under  AFOSR  Grant  94-0143  (Technical  Monitor:  J.  M.  Tishkoff)  during  the  1-year  period:  15 
February  '94-14  February  '95.  Readers  interested  in  greater  detail  than  contained  in  Section  2  are 
advised  to  consult  the  published  papers  explicitly  cited  in  Sections  2  and  5.  Copies  of  any  of  these 
published  papers  (Section  5.2  and  Appendix)  or  preprints  (Section  5.3)  can  be  obtained  by 
writing  to  the  PI:  Prof.  Daniel  E.  Rosner,  at  the  Department  of  Chemical  Engineering,  Yale 
University,  New  Haven  ,  CT  06520-8286  USA.  Comments  on,  or  examples  of,  the  applications 
of  our  research  (Section  3.4)  will  be  especially  welcome. 

An  interactive  experimental/theoretical  approach  has  been  used  to  gain  understanding  of 
performance-limiting  chemical-,  and  mass/energy  transfer-phenomena  at  or  near  interfaces.  This 
included  the  development  and  exploitation  of  seeded  laboratory  counterflow  diffusion  flame 
burners  (Section  2.1),  new  optical  diagnostic  and  particle  characterization  techniques  (Section 
2.2)  and  induction-heated  flow  reactors  (Section  2.4).  Resulting  experimental  rate  data  ,  together 
with  the  predictions  of  asymptotic  theories  ,  continue  to  be  used  as  the  basis  for  proposing  and 
verifying  simple  viewpoints  and  effective  engineering  comelations  with  a  rational  basis  for  future 
design/optimization  studies. 

2.  RESEARCH  ACCOMPLISHMENTS 

Most  of  the  results  obtained  under  Grant  AFOSR  94-0143  during  2/94-2/95  can  be 
divided  into  the  subsections  below: 

2. 1 .  FORMATION,  TRANSPORT  AND  STABILITY  OF  COMBUSTION-GENERATED  PARTICLES : 

SEEDED  LAMINAR  COUNTERFLOW  DIFFUSION  FLAME  EXPERIMENTS 

Based  on  our  recent  measurements  of  the  thermophoretic  dijfnsivity  of  flame-generated 
submicron  "soot"  particles  using  a  TiCL4(g)-seeded  low  strain-rate  countertlow  laminar  diffusion 
flame  (CDF-)  technique  (Gomez  and  Rosner,  1993)  we  showed  that  a  knowledge  of  the  relative 
positions  of  the  gas  and  particle  phase  stagnation  planes  and  the  associated  local  chemical 
environments,  can  be  used  to  control  the  composition  and  morphology  of  flame-synthesized 
particles.  These  factors  will  also  influence  particle  production  and  radiation  from  ttirbulent  non- 
premixed  "sooting"  flames. 


I 


To  obtain  fundamental  information  on  the  nucleation,  growth  and  restructuring  kinetics  of 
flame-generated  aggregates,  during  this  past  year  we  have  exploited  an  improved  "slot-type" 
burner  seeded  with  the  Ti02(s)  vapor  precursor:  titanium  tetra-iso-propoxide  (TTIP)  or  the 
Ala03(s)  vapor  precursor  trimethyl  aluminum  (TMA)  (Fig.  1)  to  carry  out  in  situ  measurements 
of  aggregate  dynamics  and  morphology  evolution.  For  the  latter  we  use  the  morphology- 
insensitive  thermophoretic  sampling  technique  (Rosner,  Mackowski  and  Garcia- Ybarra,  1991) 
with  carbon  film-coated  copper  grids  to  extract  aggregates  from  various  positions  in  the  seeded- 
CDF  for  morphological  analysis  using  transmission  electron  microscope  (TEM-)images  (Fig.  2a). 
Aggregate  data  obtained  from  CHq  flames  seeded  with  TTIP-  or  TMA-vapor  are  now  being 
obtained  and  analyzed  using  the  theoretical  methods  briefly  outlined  in  Section  2.  Apart  from 
characterizing  the  fractal  dimension,  Df,  of  aggregates  sampled  from  various  points  within  the 
CDF,  we  are  currently  examining  other  potentially  instructive  "fingerprints",  such  as  the 
distribution  of  the  angles  between  triplets  of  primary  particles,  based  on  their  projected  TEM 
images.  When  corrected  for  the  fact  that  the  centroids  of  the  particles  comprising  the  aggregate 
are  not  all  in  one  plane,  we  believe  these  histograms  will  be  useful  to  characterize  the  degree  of 
aggregate  restructuring  of  moderate  size  aggregates  in  high  temperature  flames  (see  Section  2  , 
Cohen  and  Rosner,  1993, 1994,  and  Tandon  and  Rosner,  1995).  Our  long  range  goal  is  to 
understand  how  changes  in  seed  level,  flame  stoichiometry  and  strain  rate  influence  the  nature  of 
the  particles  formed  in  this  well-defined  laminar  mixing/reaction  zone,  and  to  assess  the  relevance 
of  this  information  to  turbulent  non-premixed  reactors. 

2.2  TRANSPORT  AND  RESTRUCTURING  PROPERTIES  OF  FLAME-GENERATED 

AGGREGATES:  THEORY 

The  ability  to  reliably  predict  the  transport  properties  and  morphological  stability  of 
aggregated  flame-generated  particles  (carbonaceous  soot,  AI2O3,  Si02,...)  is  important  to  many 
technologies,  including  chemical  propulsion  and  refractory  materials  fabrication.  Needed  are 
methods  to  anticipate  coagulation  and  ultimate  deposition  rates  of  suspended  populations  of  such 
particles  in  combustion  systems.  Toward  this  end  we  are  developing  efficient  methods  for 
predicting  the  transport  properties  of  large  'fractal'  aggregates  via  a  spatially  variable  effective 
porousity  pseudo-continuum  model  .  Since  "young"  small  aggregates  are  frequently  "stringy" 
(fractal  dimensions  much  less  than  2)  we  are  also  introducing  an  aggregate  statistical 
characterization  based  on  the  pdfoi  angles  between  triplets  of  primary  particles  (Fig.  3  a,b),  and 
developing  methods  to  predict  the  evolution  of  such  pdfs,  due  to  the  restructuring  mechanisms  of 
Brownian  motion,  surface  energy  driven  viscous  flow,  and/or  capillary-condensation  (Cohen 
and  Rosner,  1993,1994).  Indeed,  these  factors  determine  the  observed  size  of  the  apparent 
"primary"  particles  comprising  soot  particles,  the  "collapse"  of  surface  area  observed  in  some 
high  temperature  systems,  and  the  associated  evolution  of  particle  transport  properties  (see,  eg., 
Tandon  and  Rosner,  1995). 

These  new  methods/results,  together  with  our  recent  estimates  of  the  spread  of  aggregate 
sizes  in  coagulating  populations,  are  now  being  used  to  predict  particle  capture  rates  on  solid 
surfaces  by  the  mechanisms  of  convective-diffusion,  turbulent  eddy-impaction,  and 
thermophoresis  (Rosner,  Tassopoulos  and  Tandon,  1993,1994).  Also  under  development  are 
methods  to  predict  interactions  between  aggregates  and  their  surrounding  vapor  environment- 
interactions  which  can  lead  to  primary  particle  growth,  or  burn-out.  Toward  this  end,  we  have 
also  developed  new  and  efficient  methods  to  predict  the  "accessible  surface  area"  of  populations 

of  aggregates  (expressed  as  a  fraction  ,  rjpop,  of  the  true  total  surface  area  (Rosner  and 
Tandon,  1994) ,  including  its  dependence  on  mean  aggregate  size  (N) ,  probing  molecule  reaction 

probability  a,  and  pressure  level  {via  Knudsen  number  based  on  primary  particle  diameter)  for 
aggregates  of  a  prescribed  structure  (fractal  dimension,  Df). 

We  are  now  carrying  out  theoretical  studies  on  the  structure  of  thin  reaction-nucleation- 
coagulation  'sublayers'  within  laminar  mixing-  and  boundary-  layers  (Fig. 7  below),  including 
stagnation  flows  similar  to  those  achieved  in  our  counterflow  burner  (Fig.  1)  and  CVD- 
impingement  flow  reactor  (  Fig.  5).  An  account  of  our  recent  studies  of  the  unusual  population 
dynamics  of  coagulating  absorbing-emitting  particles  in  strong  radiation  fields  will  be  found  in 
Aerosol  Sci.  Tech.  (Mackowski  et.ai,  1993).  For  an  overview  of  our  work  on  these  effects 
(energy  transfer  on  particle  dynamics),  see  o\xx  lEC-Research  paper:  Rosner,  et.  al,  1992. 


SPECTRUM  ANALYZER  PC 


Fig.l  "SIot"-type  laminar  counterflow  diffusion  flame  (CDF) 
burner  set-up  for  in  situ  and  extractive  experimental  studies  of 
the  nucleation,  growth,  transport  and  restructuring  of 
aggregates  in  flames  (after  Albagli,  Xing  and  Rosner,1995; 
see,  also,  Gomez  and  Rosner,  1993) 


Connection  Angle  (deg) 

Fig.  3  Distribution  of  angles  (projected)  between  touching 
'primary'  particles  in  sample  of  7  multiparticle  aggregates 
lliermophoretically  extracted  from  two  locations  on  oxidizer 
side  of  laminar  CDF  seeded  with  TTIP:  a)  5  mm  from 
flame,  b)  2  mm  from  flame  (after  Albagli,  Xing  and 
Rosner,  1994;  see,  also,  Cohen  and  Rosner,  1993,1994)) 


Fig.  2  Typical  multi-particle  aggregate  thermoplioretically 
extracted  from  laminar  CDF  (CH4/O2/N2,  equiv.  ratio: 
0.833,  nom.  strain  rate  4  s’^)  seeded  witli  Ti02  precursor 
TTIP  vapor.  TEM  image  (a)  compared  witli  ’touching  sphere’ 
idealization  (b)  (after  Albagli,  Xing  and  Rosner,  1995) 


Fig.4  Predicted  angular  dependence  of  nondimcnsional  erosion 
rate  on  an  initially  circular  metal  cylinder  in  high  Reynolds 
number  crossflow  (after  Rosner,  Tandon,  and  Labowsky 
1995).  Here  0*=66"  (incidence  angle  of  maximum  erosion 
yield),  fo=0.28  (ratio  of  erosion  yield  at  normal  incidence  to 
tliat  at  angle  0*  ,  n=2.5,  and  log-normal  population  spread, 
Og,  is  2.3) 
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Fig.  5  Axisymmelric  impinging  jet  CVD-flow  reactor  witli 
inductively  heated  "pedestal"  for  studies  of  vapor  +  particle 
co-deposition  across  laminar  BLs  (after  Rosner,  Collins  and 
Castillo,  1993,  Collins,  1994,  and  Gokoglu  et.a/.,1994) 
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Fig.  7  Predicted  local  size  distributions  of  particles  nucleated 
within  thermal  boundary  layer  near  hot  deposition  surface. 
y/5BL  =  relative  position  in  LBL;  vlv\  =  particle  volume  in 
monomer  volume  units  (after  Tandon  and  Rosner,  1995) 
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Fig. 6  Ti02(s)  deposition  rate  data  (reported  as  apparent 
heterogeneous  reaction  probability  e)  from  TTIP/O2/N2 
mainstream  showing  calculated  deposition  rate  behavior  in 
presence  of  homogeneous  reactions  to  form  non-depositing 
species  (after  Rosner,  Collins  and  Castillo,  1993,  Collins, 
1994,  Collins  and  Rosner,  1994) 
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Fig.  8  "Universar*  correlation  of  experimentally  reported  solid 
deposit  densities  based  on  a  Damkohler  number  (’’burial” 
parameter)  comparing  tlie  characteristic  times  for  surface 
reaction  and  surface  diffusion  (after  Klto,  Collins  and  Rosner, 
^994‘ 


2.3  EROSION  BEHAVIOR  OF  CERAMIC  OR  METAL  TARGETS  IN  HIGH-SPEED  ABRASIVE 

STREAMS 

By  capturing  with  simple  formulae  the  essential  features  of  available  erosion  yield 
experiments  (carried  out  using  planar  targets  with  suspended  particles  of  a  single  size),  we  have 
developed  an  efficient  method  to  predict  local  and  total  surface  erosion  rates  for  metal  or  ceramic 
targets  of  simple  geometry  exposed  to  a  'polydispersed'  population  of  abrasive  particles 
suspended  in  a  high-speed  mainstream  (Rosner,  Tandon  and  Labowsky,  1995,  Kho,  Rosner  and 
Tandon,1995,  Khalil  and  Rosner,  1995f  For  a  circular  cylinder  target  over  the  entire  Stokes 
number  range  'universal'  graphs  have  been  provided  covering  the  anticipated  dimensionless 
parameter  ranges  of  greatest  interest.  An  example  is  provided  in  Fig.  4,  which  is  a  polar  plot  of 
the  dimensionless  local  erosion  rate  of  a  metal  target  over  a  range  of  values  of  the  number-mean 
particle  size  (volume)  to  the  'critical'  size  required  to  impact  under  the  prevailing  conditions. 
Based  on  our  summary  of  available  experimental  data,  we  also  provide  representative  sets  of 
numerical  values  of  the  key  phenomenological  erosion  yield  parameters  that  emerge  from  our 
formulation,  albeit  at  Tw  =  300K.  To  deal  with  more  general,  or  'singular'  cases  of  particular 
interest  to  the  reader,  the  required  correlation  formulae  and  associated  quadrature  expressions 
have  been  provided  provided  (Rosner,  Tandon  and  Labowsky,  1995,  Kho,  Rosner  and 
Tandon,  1995,  Khalil  and  Rosner,  1995). 

Among  other  things,  the  present  model  and  our  results  make  it  clear  that  erosion  rates  will 
be  directly  proportional  to  the  abrasive  particle  mass  loading  in  the  mainstream,  and  quite  sensitive 
to  mainstream  velocity;  viz.',  even  more  sensitive  than  (where,  for  metals,  often,  n=2.5)  at 
intermediate  Stokes  numbers.  Our  present  results  also  explain  why  initially  'blunt'-nosed  metal 
objects  tend  to  become  'wedge-shaped'  (sharpened)  as  a  result  of  erosion  {cf.  Fig4  in  the  vicinity 
of  the  forward  stagnation  line),  whereas  blunt  ceramic  objects  tend  to  become  still  blunter. 


2.4  FORMATION  KINETICS  AND  MORPHOLOGY  OF  CVD-MATERIALS;  THEORY  OF  MULTI-PHASE 

BOUNDARY  LAYERS  WITH  NUCLEATION,  GROWTH  AND  THERMOPHORESIS 

A  small  impinging  jet  (stagnation  flow)  reactor  (Fig.  5)  has  been  used  to  study  the 
chemical  vapor  deposition'  (CVD-)  rates  of  refractory  layers  on  inductively  (over-)heated 
substrates  in  the  presence  of  complicating  homogeneous  reactions  of  the  vapor  precursor  (Colhns, 
1994).  These  measurements  are  being  used  to  understand  deposition  rates  and  associated  deposit 
microstructures  observed  in  highly  non-isothermal,  often  particle-containing  local  CVD 
environments.  Figure  6  shows  (logarithmic  ordinate)  our  apparent  deposition  probability  vj. 
reciprocal  surface  temperature  for  Ti02(s)  films  obtained  from  TTIP(g)  over  the  broad  surface 
temperature  range:  600-1600K.  A  mathematical  model  incorporating  finite-rate  reagent  depletion 
near  the  hot  surface  captures  the  experimentally  observed  deposition  rate  trends  (Fig.  6;  curve 

marked  Darnhom  =AhomS^/D=  8.7x104).  A  more  detailed  mathematical  model  (Tandon  and 
Rosner,  1994)  tracks  the  size  distributions  of  nucleated  particles  within  such  BLs  (Fig.  7), 
including  the  particle-i-vapor  "co-deposition"  rate  at  the  surface  y=0,  and  the  associated  surface 

roughness.  ... 

Missing  from  previous  predictions  of  film  deposition  rates  has  been  essential  information 
about  deposit  microstructure  and  ancillary  thermophysical  properties.  Toward  this  end  we  have 
investigated  the  possibility  of  a  rational  correlation  of  vapor  deposit  grain  densities  based  on  the 
notion  that  close-packed  dense  deposits  should  be  possible  only  when  the  characteristic  time  for 
adatom(molecule)  surface  diffusion  is  short  compared  to  the  time  for  surface  reaction  and 
"burial".  Our  preliminary  results  (Fig.  8;  after  Kho,  Collins  and  Rosner,1994)  are  encouraging 
and  the  addition  of  further  relevant  characteristic  times  to  our  scheme  may  open  the  door  to  the 
corelation/'anticipation’  of  other  important  deposit  characteristics. 


In  our  OSR-sponsored  Yale  HTCRE  Lab  research  during  this  program,  briefly  reviewed 
here,  we  have  shown  that  new  methods  for  rapidly  measuring  particle  transport  properties  and 
rates,  combined  with  recent  advances  in  boundary  layer  theory,  provide  useful  means  to  identify 
and  incorporate  important,  but  often  previously  neglected,  mass  transport  phenomena  in  many 
multiphase  propulsion  engineering  and  materials  engineering  design/optimization  calculations. 

Despite  formidable  complexities  to  be  overcome  in  the  design  and  operation  of  air- 
breathing  propulsion  power  plants  utilizing  a  broad  spectrum  energetic  fuels  these  particular 
techniques  and  results  are  indicative  of  the  potentially  useful  simplifications  and  generalizations 
which  have  emerged  from  this  program's  fundamental  AFOSR-funded  research  studies  of 
combustion-generated  particle  transport  mechanisms,  deposition  ^d/or  erosion.  It  is  hoped  that 
this  Annual  Report  and  its  supporting  (cited)  papers  will  facilitate  the  refinement  and/or 
incorporation  of  some  of  the  present  ideas  into  engineering  design  procedures  of  much  greater 
generality  and  reliability.  This  work  has  already  helped  identify  new  directions  where  research 
results  could  have  a  significant  impact  on  engineering  practice  in  both  the  defense  and  civilian 
sectors  of  the  US  economy  (Section  3.4). 


3.  ADMINISTRATIVE  INFORMATION: 

PERSONNEL,  PRESENTATIONS,  APPLICATIONS, 

"COUPLING’  ACTIVITIES 

The  following  sections  summarize  some  pertinent  'non-technical'  facets  of  the 
abovementioned  Yale  HTCRE  Lab/AFOSR  research  program  during  the  period;  2/15/94-2/14/95: 

3.1  Personnel  j 

The  present  results  (Sections  2  and  5)  are  due  to  the  contributions  of  the  individuals  hsted 
in  Table  3.1-1,  which  also  indicates  the  role  of  each  researcher  and  the  relevant  time  interval  of 
the  activity.  It  will  be  noted  that,  in  addition  to  the  results  themselves,  this  program  has 
simultaneously  contributed  to  the  research  training  of  a  number  of  students  and  2  recent  PhDs, 
who  will  now  be  in  an  excellent  position  to  make  future  contributions  to  technologies  oriented 
toward  air-breathing  chemical  propulsion,  and  high-tech  materials  processing. 


Table  3.1-1  Summary  of  Research  Participants^  on  AFOSR  Grant :  94-0143 

TRANSPORT  PHENOMENA  AND  INTERFACIAL  KINETICS 
IN  MULTIPHASE  COMBUSTION  SYSTEMS 


Name 

Status^ 

Principal  Research  Activitvh 

Albagli,  D. 

PDRA 

particle  prod/char,  in  CDFs 

Castillo,  J.L. 

VS 

reactions  in  thermal  BLs 

Collins,  J. 

GRA 

CVD  of  ceramic  coatings 

Garcia- Ybarra,  P. 

VS 

aggregate  transport  theory 

Khalil,  Y.F. 

RA 

erosion  of  ceramics 

Kho,  T. 

GRA 

correl.  of  density  of  CVD  coatings 

Konstandopoulos,  A.G. 

VS 

combined  inertia  +  thermophoresis 

Koylu,  U. 

PDRA 

particle  prod/char,  in  CDFs 

Labowsky,  M.  J. 

VS 

inertial  impaction  and  erosion  of  metals 

Papadopoulos,  D 

GRA 

boundary  conditions  at  G/S  interfaces 

Rosner,  D.E. 

PI 

program  direction-dep.  theory/exp 

Tandon,  P. 

GRA 

transport  phenom.  in  BLs  and  CDFs 

Xing,  Y. 

GRA 

particle  prod/char.exps.  in  CDFs 

^  PDRA=Posi-doctoral  Research  Asst  GRA=  Graduate  Research  Assistant 
PI  =  Principal  Investigator  VS  =  Visiting  Scholar  RA  =  Research  Affiliate 

See  Section  5  for  specific  references  cited  in  text  (Section  2) 
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3.2  Cooperation  with  US  Industry 

The  research  summarized  here  was  supported  by  AFOSR  under  Grant  94-0143  (2/15/94- 
2/14/95).  In  these  fields,  the  Yale  HTCRE  Laboratory  has  also  been  the  beneficiary,  over  the 
past  decade,  of  smaller  grants  from  many  interested  U.S.  industrial  corporations,  including 
groups  within  GE,  DuPont,  Union  Carbide  (now  Advanced  Ceramics  Corp.)  and  Shell,  as  well 
as  the  feedback  and  advice  of  principal  scientists/engineers  from  each  of  these  corporations  and 
Combustion  Engineering- ABB  and  Textron.  We  appreciate  this  level  of  collaboration,  and  expect 
that  it  will  accelerate  inevitable  applications  of  our  results  in  areas  relevant  to  their  technological 
objectives  (see,  also.  Section  3.4,  below).  We  also  expect  these  interactions  to  intensify  as  a 
result  of  the  recent  formation  of  an  interdisciplinary  Combustion  Research  Center  at  Yale. 

3.3  Presentations  and  Research  Training 

Apart  from  the  publications  itemized  in  Section  5  and  our  verbal  presentation  (of  progress) 
at  the  AFOSR  Contractors  Meeting  (June  8-10,  Tahoe,  NV) ,  our  results  over  this  past  year  have 
also  been  presented  at  some  10  seminars/conferences-including  annual  or  topical  conferences  of 
the  following  professional  organizations; 

Particle  Technology  Forum  (8/94,  Denver  CO)) 

AIChE  (11/94,  San  Francisco  CA) 

4th  Int.  Aerosol  Conference  (Aug.29-Sept.2,1994  @  UCLA,  CA) 

In  addition,  during  the  period:  2/15/94-2/14/95,  the  PI  presented  seminars  at  the  following 
Laboratories:  ALCOA,  PPG,  and  Technical  Universities:  Istanbul,  Athens,  Thessaloniki,  ETH- 
Zurich  ,  Budapest  (these  were  a  portion  of  Prof.  Rosner's  Fall  1995  Yale  Academic  Leave 
itinerary). 

This  program  involved  completion  of  the  PhD  dissertation  research  of  two  Yale  ChE 
Dept,  graduate  students  (Joshua  Collins,  who  defended  on  3/4/94,  and  Pushkar  Tandon,  who 
defended  on  2/3/95  and  will  remain  at  Yale  beyond  2/15/95  as  a  Post-doctoral  Research 
Engineering  Scientist  (cf.  Table  3.1-1). 


3.4  Some  Known  Applications  of  Recent  Yale-HTCRE  Lab  Research  Results 

It  has  been  particularly  gratifying  to  see  direct  applications  of  some  of  this  generic 
AFOSR-supported  particle  and  vapor  mass  transfer  research  in  more  applications-oriented 
investigations  reported  in  recent  years.  Indeed,  the  writer  would  appreciate  it  if  further  examples 
known  to  the  reader  can  be  brought  to  his  attention.  Not  surprisingly,  because  of  the  inevitable 
"induction  period"  following  publication,  our  very  recent  erosion  rate  studies  will  be  covered  in 
next  year's  report. 

Our  AFOSR  supported  research  on  soot  deposition  rates  from  flowing  laminar  or 
turbulent  combustion  gases  has  been  applied  by  Aerojet  Corp.  (D.  Makel  et.al.  A990)  to  develop 
a  model  for  application  to  rocket  chambers  and  nozzles  (with  NASA  support).  Extensions  to  jet 
engine  nozzles  are  currently  being  made  by  M.T.  Nys  at  Pratt  and  Whitney  Engine  Business  in 
W.  Palm  Beach  FL. 

In  the  area  of  multicomponent  vapor  deposition  in  combustion  systems  applications  of  our 
predictive  methods  (for  "chemically  frozen"  (Rosner  et.al,  1979)  and  LTCE  multicomponent 
laminar  boundary  layers)  have  been  made  by  British  Coal  Corporation-Power  Generation  Branch 
(I.  Fantom,  contact)  in  connection  with  their  topping  cycles  which  run  gas  turbines  on  the 
products  of  fluidized  bed  coal  combustors/gasifiers.  Also,  in  combustion  research  many  groups 
(eg.,  Dobbins  et.al.  (Brown  U.),  Faeth  et.al.  (U.  Mich.),  Katz  et  al.  (J.  Hopkins  U.))  are  now 
utilizing  "thermophoretic  sampling"  techniques  to  exploit  the  size-  and  morphology-insensitive 
capture  efficiency  characteristics  that  we  have  proven  in  our  AFOSR  research  (Section  2.1). 

Our  AFOSR  and  NASA  fundamental  research  on  chemical  element  segregation  in  the 
CVD  of  refractory  ceramics  (eg.,  SiC  and  metal  borides)  (see,  eg.,  Collins  and  Rosner,  1991, 
1992)  is  evidently  of  use  to  AFML  contractors  synthesizing  controlled  stoichiometry  fibers  for 
light  weight/high  strength  composites  (Americom,  Textron). 
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For  calculating  suspended  particle  concentrations  along  trajectories  outside  of  aircraft 
(involved  in  atmospheric  sampling),  or  inside  of  CVD  reactors,  A.  S.  Geller  and  D.  J.  Rader  of 
Sandia- Albuquerque  have  adopted  a  method  developed  in  our  earlier  AFOSR  work  (Fernandez 
de  la  Mora,  1981),  and  recently  applied  in  our  own  studies  of  particle  motion  in  laminar  boundary 
layers  with  streamwise  curvature  (Konstandopoulos  and  Rosner,  1995). 

Ongoing  work  at  MIT  (Walsh  et.al.  1992),  PSI  and  Sandia  CRF  has  incorporated  our 
rational  correlation  of  inertial  particle  impaction  {e.g.  a  cylinder  in  cross-flow)  in  terms  of  our 
effective  Stokes  number  (Israel  and  Rosner,  1983,  and  Konstandopoulos  et.  al.  1993).  Recent 
applications  of  our  AFOSR  and  DOE-supported  research  (on  the  correlation  of  inertial  impaction 
by  cylinders  in  crossflow)  have  also  been  made  by  the  National  Engineering  Laboratory  (NEL)  of 
Glasgow  Scotland  .  NEL  is  apparently  developing  mass-transfer  prediction  methods  applicable  to 
waste-heat  recovery  systems  in  incinerators,  as  well  as  pulverized  coal-fired  boilers.  These 
applications  are  somewhat  similar  to  those  reported  by  the  Combustion  R&D  group  at  MIT  and 
Penn  State  U,  and  are  now  being  taken  up  by  VTT-Energy/ Aerosol  Technology  Group,  in 
Finland. 

Explicit  use  of  our  studies  of  self-regulated  "capture"  of  incident  impacting  particles 
(Rosner  and  Nagaragan,  1987)  is  being  made  in  current  work  on  impact  separators  and  ceramic 
heat  exchangers  for  coal-fired  turbine  systems  in  high  performance  stationary  power  plants.  Other 
potential  applications  arise  in  connection  with  "candle  filters"  used  to  remove  fines  (sorbent 
particles,...)  upstream  of  the  turbines.  A  useful  summary  of  work  in  these  interrelated  areas 
(Solar  Turbines,  Textron  Defense  Systems,  Hague  International,...)  was  presented  at  the 
Engineering  Foundation  Conference  Inorganic  Transformations  and  Ash  Deposition  During 
Combustion.,  the  proceedings  of  which  appeared  in  1992. 

Clearly,  fruitful  opportunities  for  the  application  of  our  recent  "non-Brownian"  convective 
mass  transfer  research  now  exist  in  many  of  the  programs  currently  supported  by  the  US  Air 
Force,  as  well  as  civilian  sector  R&D.  Indeed,  based  on  the  abovementioned  10  papers  we  have 
prepared  (now  in  press  or  submitted  (Section  5.3))  as  well  as  our  lecture-presentations  (Section 
3.3) ,  we  expect  that  these  applications  will  multiply  rapidly. 

4.  CONCLUSIONS 

In  the  OSR-sponsored  Yale  HTCRE  Lab  research  during  the  period:  2/15/94-2/14/95, 
briefly  described  above,  we  have  shown  that  new  methods  for  rapidly  measuring  particle-mass 
transport  properties  and  rates,  combined  with  our  recent  advances  in  mass  transport  theory, 
provide  useful  means  to  identify  and  incorporate  important,  but  previously  neglected,  mass 
transport  phenomena  in  many  chemical  propulsion  engineering  and  materials  engineering 
design/optimization  calculations.  Among  our  research  results  described  in  detail  in  the  cited 
references  (Section  5),  perhaps  the  most  noteworthy  are  the  development/reporting  of: 

R1  rational  methods  to  predict  the  restructuring  kinetics,  and  transport  property  evolution 
for  aggregated  'soot'  particles  in  high  pressure  combustion  gases 
R2  rational  methods  to  predict  the  erosion  rates  of  metal  and  ceramic  solid  surfaces  in 
abrasive,  high-speed  streams 

R3  Experimental  data  and  quantitative  methods  for  predicting/correlating  the  effects  particle 
formation  in  the  CVD  thermal  boundary  layei/on  the  rate  and  quality  of  vapor- 
deposited  ceramic  thin  films  ' 

R4  Damkohler  number-based  rational  correlation  for  the  density  of  vapor  deposited  solid 

Despite  formidable  complexities  to  be  overcome  in  the  design  and  operation  of  mobile 
(and  stationary)  power  plants  utilizing  a  broad  spectrum  of  energetic  fuels  the  abovementioned 
techniques  and  results  (Sections  2,5)  are  indicative  of  the  potentially  useful  simplifications  and 
generalizations  emerging  from  our  present  fundamental  AFOSR-funded  research  studies  of 
combustion-generated  particle  transport  mechanisms  and  interfacial  reactions  relevant  to  the 
synthesis  of  refractory  materials.  It  is  hoped  that  this  Annual  Report  and  its  supporting  papers 
(Section  5)  will  facilitate  the  incorporation  of  many  of  the  present  ideas  into  design  and  test 
procedures  of  greater  generality  and  reliability.  This  work  has  also  helped  identify  new  directions 
where  it  is  anticipated  that  research  results  from  this  AFOSR  program  have  a  significant  impact  on 
future  USA  DOD  and  civilian  sector  engineering  practice. 
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ABSTRACT 

We  model  a  cold-wall  atmospheric  pressure  impinging  jet  reactor 
to  study  the  CVD  and  gas-phase  nucleation  of  TiOj  from  a  titanium 
tetra-iso-propoxide  (TTIP) /oxygen  dilute  source  gas  mixture  in 
nitrogen.  The  mathematical  model  uses  the  computational  code  FIDAP 
and  complements  our  recent  asymptotic  theory  for  high  activation 
energy  gas-phase  reactions  in  thin  chemically  reacting  sublayers. 
The  numerical  predictions  highlight  deviations  from  ideality  in 
various  regions  inside  the  experimental  reactor.  Model  predictions 
of  deposition  rates  and  the  onset  of  gas-phase  nucleation  compare 
favorably  with  experiments.  Although  variable  property  effects  on 
deposition  rates  are  not  significant  (-11%  at  lOOOK) ,  the  reduction 
of  rates  due  to  Soret  transport  is  substantial  (-75%  at  lOOOK) . 


1.  INTRODUCTION 

Production  of  sophisticated  materials  with  superior  properties 
by  CVD  requires  an  understanding  of  coupled  transport  and  chemical 
processes.  This  goal  can  be  realized  in  research  programs  by 
combining  both  experiments  and  modeling. 

Recent  Yale  University  research  focused  on  TiOj  CVD  with  simul¬ 
taneous  gas-phase  reaction  leading  to  particle  nucleation  [1-2]. 
A  simple  asymptotic  theory  treated  the  onset  of  reactions  leading 
to  particle  nucleation  and  reduced  CVD  rates  at  high  surface  tempe¬ 
ratures.  The  theory  assumed  that  all  gas-phase  reactions  are  rest¬ 
ricted  to  a  thin  chemical  sublayer  adjacent  to  a  hot  CVD  surface 
due  to  their  large  activation  energies.  They  also  developed  a  lab- 
scale,  cold-wall,  impinging  jet  CVD  reactor  [3].  Using  experiments 
and  appropriate  asymptotic  and  numerical  models,  we  expect  to  guide 
future  reactor  designs  and  scale-up,  resulting  in  maximum  CVD  rates 
while  avoiding  harmful  particle  nucleation. 

The  present  axisymmetric ,  numerical  model  points  to  the  role  of 
bouyancy  and  recirculation  in  the  experimental  reactor  and  clari¬ 
fies  the  effect  of  Soret  transport  in  modifying  TTIP  mass  transfer 
rates.  It  provides  the  ability  to  capture  the  transition  from 
surface  kinetics  to  gas-phase  transport  control.  Furthermore,  it 
extends  the  earlier  asymptotic  theory  [1]  by  removing  the  restric¬ 
tion  that  the  activation  energies  of  gas-phase  reactions  should  be 
"large”  (i.e.  the  chemical  sublayer  should  be  "thin") . 

For  the  case  studied  at  Yale  [1-3]  for  the  CVD  of  Ti02  from  di¬ 
lute  TTIP/Oj  in  Nj,  the  ratio  of  chemical  sublayer  to  thermal  bound¬ 
ary  layer  thickness  is  estimated  to  be  1/5  at  the  onset  of  particle 
nucleation  (T=^1050K)  ,  which  may  be  too  large  for  the  accuracy  of 
the  asymptotic  theory.  Hence,  our  numerical  model  may  ultimately 
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be  better  able  to  capture  the  effects  of  homogeneous  chemistry 
leading  to  an  observed  deposition  rate  fall-off  above  -1050K. 


2.  EXPERIMENT 

We  use  a  cold-wall,  atmospheric  pressure,  axisymmetric,  imping- 
ing-jet  reactor  (Fig.  1).  The  liquid  TTIP  source  is  a  constant 
temperature  bubbler.  After  mixing  TTIP  vapor  with  excess  oxygen 
(to  "burn  away"  any  co-deposited  carbon)  and  diluting  with  nitrogen 
in  a  short  mixing  chamber,  the  gas  jet  emerges  from  a  converging 
nozzle  and  impinges  on  a  polished  quartz  substrate  (diam.  1.3cm). 
The  substrate  and  the  alumina  substrate  holder  are  supported  from 
below  by  an  RF-heated  graphite  susceptor.  The  susceptor  tempera¬ 
ture  is  measured  using  a  Pt-30%Rh/Pt-6%Rh  thermocouple  just  below 
the  substrate.  This  thermocouple  reading  is  correlated  with  direct 
substrate  surface  temperature  measurements  made  in  situ  using  known 
melting  point  lacquers  (with  -3-4%  accuracy) .  Gas  and  surface  tem¬ 
peratures  at  several  other  locations  are  also  measured  by  thermo¬ 
couples.  The  concentration  of  water  vapor  in  the  reactor  is  less 
than  3ppm.  Deposition  rates  are  measured  by  in  situ  interfero¬ 
metry  and  confirmed  by  ex  situ  weight  gain.  Further  details  of 
the  experimental  system  and  operating  procedure  are  given  in  [3]. 


3.  PHYSICO-CHEMICAL/NUMERICAL  MODEL 

Our  modeling  study  adopts  the  finite-element-based  computational 
fluid  dynamics  code  FIDAP.  The  new  version  of  FIDAP  incorporates 
many  phenomena  relevant  to  CVD  processes,  such  as  temperature 
dependent  fluid  density,  transport  and  thermodynamic  properties, 
Soret  diffusion,  and  gas-phase  and  surface  chemical  reactions. 

We  approximate  the  reactor  by  using  a  2-D  axisymmetric  geometry 
(Fig.  2).  The  inlet  and  exit  are  placed  sufficiently  far  from  the 
deposition  zone  to  eliminate  numerical  uncertainties  on  rate 
predictions.  The  numerical  mesh  above  the  substrate  is  fine  enough 
to  resolve  the  chemical  sublayer.  Temperature  boundary  conditions 
are  interpolations  of  measured  surface  temperatures  in  the  Yale 
axisymmetric  impinging  jet  reactor. 

We  assume  that  the  gas-phase  decomposition  of  TTIP  is  first  or¬ 
der  with  a  rate  constant  of  3 . 96xl0^exp  (-8480/T)  s'^  (T  is  in  Kel¬ 
vins)  [4].  We  also  assume  that  the  gas-phase  products  form  nonde¬ 
positing  particles.  We  fit  our  low  substrate  temperature  rate  data 
to  a  Soret-corrected  Arrhenius  expression,  1. 21xlo''exp(-16480/'r) 
m/s,  to  obtain  a  pseudo-first  order  surface  reaction  rate  constant. 
This  expression  corresponds  to  a  TTIP  reactive  sticking  probability 
of  unity  at  987K.  For  simplicity,  the  surface  reaction  rate 
constant  is  kept  at  its  987K  value  for  higher  substrate  tempera¬ 
tures  so  that  the  sticking  probability  does  not  exceed  unity. 

The  Soret  diffusion  factor  for  dilute  TTIP  in  Nj  is  calculated 
from  kinetic  theory  (the  estimated  Lennard-Jones  parameters  for 
TTIP  are  a=8.13A  and  e/k=589K)  and  fitted  to  the  expression 
aT=1.971[l-(223.6)/T]  .  Unfortunately,  FIDAP  is  currently  restricted 
to  constant  Soret  dif fusivities,  Dy=pDwaT,  hence  we  evaluate  the  gas 
(N2)  density  p,  TTIP  Fickian  diffusivity  D,  and  Oj  at 
T,i,„=(T„,,„e+T,,rtace)/2,  and  USB  the  TTIP  mass  fraction  a.  at  the  T,„, 
location  for  the  corresponding  diffusion  limit  calculation. 


172 


4.  RESULTS  AND  DISCUSSION 


First,  we  compare  FIDAP  predictions  to  the  observed  flow  field 
inside  the  reactor.  The  experiments  use  argon  seeded  with  fine  ti- 
tania  particles  for  flow  visualization.  Figure  3  shows  a  typical 
comparison  at  a  Reynolds  number  of  750  based  on  nozzle  diameter  at 
STP  and  900K  substrate  temperature.  Indeed,  FIDAP  can  demonstrate 
recirculations  a)  as  the  jet  emerges  from  the  nozzle  due  to  flow 
separation,  and  b)  on  the  sides  of  the  hot  susceptor  due  to  bouyan- 
cy-induced  convection.  These  recirculations  are  not  detrimental 
because  from  the  photograph  one  can  see  that  particles  trapped  in 
the  recirculation  regions  do  not  diffuse  into  the  reagent  jet,  and 
we  estimate  that  reaction  product  vapor  species  are  also  unable  to 
penetrate  appreciably  into  the  jet.  The  Richardson  number  based 
on  nozzle  to  substrate  distance  is  -0.07  for  this  case.  Lower  flow 
rate  (Re<100)  measurements  at  higher  substrate  temperatures  are 
hampered  by  buoyancy  effects  in  the  jet  above  the  substrate.  We 
discuss  elsewhere  [5]  the  operating  conditions  needed  to  avoid 
bouyancy-induced  convection  in  such  reactors.  Our  numerical 
analysis  is  capable  of  handling  such  non-ideal  behavior. 

For  the  results  discussed  below  we  fix  the  gas  flow  rate  at  5160 
seem  (Re“500  for  Nj)  and  the  substrate-temperature-based  TTIP 
concentration  at  05=2.5x10'®  Kgmol/m\  Deposition  rates  are  repeated 
as  effective  reaction  probabilities  e  defined  as  e=n"/ [ ( 1/4 ) vC^]  , 
where  n"  and  v  are  the  molar  flux  and  mean  thermal  speed  of  TTIP, 
respectively,  evaluated  at  the  surface. 

Figure  4  depicts  Soret  diffusion  and  temperature-dependent  gas 
transport-property  effects  on  predicted  rates  at  the  stagnation 
point  for  transport  controlled  conditions.  The  rate  reduction  due 
to  Soret  transport  is  >75%,  whereas  the  effect  of  variable 
properties  is  <15%,  for  T>1000K. 

The  model  predicts  reasonable  deposition  rates  over  the  entire 
regime  from  surface  kinetics  to  transport  control  by  incorporating 
the  surface  kinetics  extracted  from  our  experiments  (Fig.  5) .  With 
the  gas-phase  kinetics  of  [4],  the  model  agreement  with  experiment 
is  less  satisfactory,  though  with  similar  trends.  As  the  depletion 
of  TTIP  by  the  gas-phase  reaction  increases  (i.e.  chemical  sublayer 
gets  thicker)  at  higher  substrate  temperatures,  thereby  decreasing 
the  TTIP  surface  flux  (CVD  rate),  FIDAP's  handling  of  Soret  trans¬ 
port  becomes  more  inconsistent.  This  is  because  the  constant  D^, 
used  by  FIDAP  as  defined  above  for  the  calculation  of  Soret  trans¬ 
port,  does  not  exhibit  the  expected  diminishing  Soret  contribution 
as  TTIP  concentration  diminishes.  Hence,  our  rate  predictions  with 
Soret  effect  above  substrate  temperatures  of  1200K  are  not  shown. 
The  cause (s)  of  our  underprediction  of  the  steepness  of  the  high 
temperature  rate  fall-off  is  currently  under  investigation. 

The  onset  of  homogeneous  nucleation  is  inferred  from  a  sudden 
drop  in  deposition  rate  at  a  certain  substrate  temperature.  The 
reliability  of  the  asymptotic  theory  prediction  of  this  temperature 
depends  on  whether  the  chemical  sublayer  is  indeed  "much"  thinner 
than  the  thermal  boundary  layer  at  the  prevailing  conditions. 
Figure  6  depicts  the  TTIP  mass  fraction  profiles  calculated  by 
FIDAP  at  the  stagnation  point  at  a  substrate  temperature  of  1200K 
with  and  without  gas-phase  reaction.  We  also  plot  the  difference 
between  the  two  mass  fraction  curves  to  demonstrate  the  extent  of 
gas-phase  reaction,  as  well  as  the  corresponding  temperature  pro¬ 
file.  The  predicted  rate  at  this  temperature  with  gas-phase  che- 
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FIGURE  6:  The  effect  of  gas  phase  reaction  on  mass  fraction  FIGURE  7:  Radial  variation  of  predicted  deposition  rates  across 

profile  at  1200K.  Soret  effect  is  included.  susceptor  surface. 


is  obsGirvalDly  lowsir  thcin  ths  on©  without  g3s— phssG  chGinist— 
ry.  If  the  peak  of  this  curve  is  taken  as  a  measure  of  the  thick¬ 
ness  of  the  chemical  sublayer,  it  is  substantial  (30%)  compared 
with  the  thermal  boundary  layer  thickness.  The  amplitude  and  width 
of  the  difference  curve  would  be  larger  at  even  higher  substrate 
temperatures.  Therefore,  under  such  conditions  predictions  of  the 
full  numerical  approach  are  expected  to  be  more  reliable  than 
predictions  of  the  simple  asymptotic  theory.  Also  noteworthy  is 
the  slight  "enrichment”  in  the  mass  fraction  profiles  away  from  the 
surface  due  to  the  Soret  effect. 

Figure  7  shows  that  the  predicted  rates  on  the  susceptor  are  not 
uniform  in  the  radial  direction  beyond  the  substrate  radius  (6.35mm 
vs.  16mm  for  the  susceptor) .  However,  because  the  deposition  rate 
is  roughly  uniform  across  the  substrate,  we  have  compared  rate 
measurements  to  stagnation  point  rate  predictions  (Figs.  4  &  5)  . 


5.  CONCLUSIONS 

Our  numerical  model  can  correctly  describe  the  flow  field  in  our 
experimental  atmospheric  pressure,  impinging  jet  reactor  used  for 
studying  the  CVD  and  gas-phase  nucleation  of  titania.  The  model 
points  to  some  buoyancy  and  flow  recirculation  effecpts  inside  the 
hot  substrate  reactor  affecting  the  "ideal"  stagnation  point  flow 
behavior.  It  can  guide  the  interpretation  of  deposition  rate  mea¬ 
surements  over  a  substrate  temperature  range  (600-1400K)  covering 
surface-kinetics-,  gas-phase  transport-,  and  gas-phase  reaction- 
governed  regimes.  Soret  transport  effects  on  predicted  rates  [6] 
are  significant  (-75%  reduction  at  lOOOK) ,  but  variable  property 
effects  are  modest  (<15%) .  Our  numerical  model  extends  the  capabi¬ 
lity  of  the  earlier  asymptotic  theory  to  predict  the  onset  temper¬ 
ature  of  nucleation  and  the  associated  reduced  deposition  rates  to 
include  conditions  where  the  chemical  sublayer  is  not  "much"  thin¬ 
ner  than  the  thermal  boundary  layer.  Future  work  will  include 
effects  of  flow  rate  and  TTIP  concentration  [7],  and  a  more  careful 
examination  of  the  high  temperature  rate  fall-off  regime. 
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Aggregates  (composed  of  targe  numbers  of  primary  particles)  are  produced  in 
many  engineering  environments.  One  convenient  characterization  is  the  fractal  dh 
mension,  the  exponent  describing  how  the  number  of  primary  particles  in  each 
aggregate  scales  with  radial  distance  from  its  center  of  mass.  We  describe  a  finite- 
analytic  pseudo-continuum  prediction  of  the  normalized  accessible  surface  area  of 
on  isothermal  quasi-spherical  fractal  aggregate  containing  N  (»  J)  primary  par¬ 
ticles,  on  the  surfaces  of  which  a  first-order  chemical  process  occurs.  Results  are 
displayed  for  specific  fractal  dimensions  (2.5,  2.18,  and  1 .8)  frequently  observed 
in  aggregating  systems.  An  effective  Thiele  modulus  is  used  to  develop  an  efficient 
and  accurate  scheme  for  predicting/correlating  the  effectiveness  factor  for  an  ag¬ 
gregate  containing  N  primary  particles  in  terms  of  aggregate  fractal  dimension, 
reaction  probability,  and  Knudsen  number.  Our  methods  now  allow  calculations  of 
the  accessible  surface  area  of  populations  of  aggregates,  provided  pdf-{N,  Dj, 
..  .y  is  known  for  the  populations  of  interest. 


Introduction 

In  many  technologies  involving  suspensions  the  dispersed 
phase  is  comprised  of  aggregated  or  agglomerated  particles, 
themselves  comprised  of  a  large  number,  N,  of  primary  par¬ 
ticles  or  spherules  which  have  coagulated  (usually  as  the  result 
of  their  Brownian  motion).  While  the  resulting  aggregates  may 
be  quasi-spherical  and  compact,  with  relatively  sharp,  if  ir¬ 
regular  boundaries,  more  often  they  are  sparse,  apparently 
random  clusters  of  spherules  quite  permeable  to  the  surround¬ 
ing  solvent  and  its  molecular  solutes  (Figure  1).  Probably  the 
most  studied  example  is  that  of  carbonaceous  soot  formed 
deliberately  or  inadvertently  in  the  combustion  of  hydrocarbon 
fuels  using  air  as  a  gaseous  oxidizer  (Medalia  and  Hickman, 
1960;  Siegla  and  Smith,  1981).  How'ever,  qualitatively  similar 
aggregated  smokes  arise  in  the  vapor  phase  synthesis  of  in¬ 
organic  pigments  [for  example,  TiOjfs)]  or  viscosity  modifiers 
(for  example,  fumed  silica)  (Ulrich,  1984;  Megaridis  and  Dob¬ 
bins,  1990),  and  in  the  combustion  of  metals  (for  example,  in 
flares  or  chemical  flash  lamps),  or  the  production  of  catalytic 
materials  (Formenti  ei  al.,  1972). 

Predictions  of  the  performance  of  systems  containing  such 
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aggregates  require  a  quantitative  understanding  of  the  laws 
governing  the  exchange  of  mass-,  energy-,  and/or  momentum 
between  each  aggregate  and  the  prevailing  host  fluid.  In  par¬ 
ticular,  we  address  here  the  question  of  w'hai  fraction  of  the 
surface  area  of  the  primary  spherules  is  in  fact  accessible  to  a 
diffusing  solute  molecule  in  the  “continuous’*  solvent  phase— 
a  high  temperature  combustion  gas  mixture  in  the  above  men¬ 
tioned  examples.  For  example,  in  the  case  of  carbonaceous 
soot  aggregates  this  fundamental  question  arises  in,  say,  cal¬ 
culating  the  growth  rate  of  the  primary  panicles  due  to  C2H2(g) 
molecule  impacts,  the  adsorption  rate  of  possibly  carcinogenic 
condensible  molecules,  or  the  ultimate  oxidation  rate  of  the 
aggregates  due  to  such  combustion  gas  species  as  OH(g),  0(g), 
02(g),  C02(g),  and  so  on.  For  more  volatile  materials  the 
present  theory  could  also  be  used  to  estimate  large  aggregate 
sublimation  rates. 

For  5 AT? ^7// aggregates,  comprised  of  only  a  few  primary  par¬ 
ticles  (say  2<jV<20)  in  a  specified  geometrical  configuration 
(for  example,  linear  chain,  or  quasi-spherical  packing  of  uni¬ 
form  density,  and  so  on)  it  is  no  longer  impractical  to  make 
detailed  calculations  of  transport  phenomena  in  the  complex 
shaped  space  exterior  to  the  primary  panicles  (Rosner  et  al., 
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Figure  1.  Simulated  cluster-monomer  ballistically  grown 
three-dimensional  aggregate  containing  10^ 
primary  particles. 

(After  P.  Meakin,  as  reported  in  Schaefer  cl  al.,  1988)  {D,^  2  5). 


Figure  2.  Mathematical  model  of  a  quasi  spherica!  ag¬ 
gregate  comprised  of  W-primary  spherules 
showing  the  important  characteristic  radii:  R, 
(primary  particle),  (aggregate),  and  R^,, 
(characteristic  aggregate  outer  radius  (Eq.  4)), 


1991;  Mackowski,  1990,  1994),  however,  this  approach  be¬ 
comes  tedious  and  error-prone  for  random  aggregates  com¬ 
prised  of  rather  more  particles,  perhaps  10^<A^<  10'',  as  have 
been  obseived  in  coagulation-aged  It  is  precisely  for 

populations  containing  very  large  aggregates  in  high-pressure 
gases  or  liquids  that  the  above  mentioned  shielding  effects  of 
exterior  particles  is  expected  to  become  appreciable,  reducing 
their  accessible  surface  area  well  below  that  of  the  sum  of  all 
the  primary  particles  present. 

By  viewing  each  large  quasi-spherical  aggregate  as  a  porous 
reactive  object  with  a  suitable  radial  variation  of  effective 
(pseudo-homogeneous)  rate  constant  and  effective  reagent  dif¬ 
fusion  coefficient  we  seek  and  provide  convenient  results  for 
the  accessible  area  of  large  aggregates,  suitable  for  future  cal¬ 
culations  of  mass  transfer  between  populations  of  such  ag¬ 
gregates  and  the  coexisting  carrier  fluid,  here  assumed  to  be 
a  vapor  mixture  with  non-negligible  molecular  mean-free-path. 
The  general  problem  is  showm  to  be  conceptually  identical  to 
that  of  predicting  the  effectiveness  factor  of  much  larger  solid 
catalyst  (support)  pellets  in  a  fixed-bed  chemical  reactor  (Aris, 
1975;  Butt,  1980;  Satterfield,  1970;  Rosner,  1986,  1990)  except 
that  we  must  now'  also  embrace  the  possibility  of  quasi-spher- 
ical  aggregates  that  are  so  sparse  as  to  be  essentially  transparent 
to  probing  molecules  traveling  on  straight  line  (balllistic)  tra¬ 
jectories.  Each  class  of  possibilities  is  conveniently  character¬ 
ized  by  the  so-called  fractal  dimension  or  scaling  exponent, 
Dj,  defined  by  d In  N/d In  r,  where  N  is  the  number  of  primary 
particles  contained  in  an  imaginary  sphere  of  radius  r  drawn 
about  the  aggregate  center  of  mass  (Figure  2)  (Meakin,  1983; 
Mountain  et  al.,  1984;  Dobbins  et  al.,  1987;  Schmidi-Ott  et 
al.,  1988a,  1990).  In  the  absence  of  aggregate  restructuring, 


Dj  has  been  observed  to  be  about  1.8  in  a  wide  variety  of 
systems  (physically  or  optically  sampled,  as  well  as  computer- 
simulations  of  aggregate-aggregate  encounters),  but  probably 
because  of  partial  restructuring,  aggregates  have  been  reported 
with  intermediate  Dy-values,  for  example,  2.18  for  (he  silver 
aggregates  experimentally  studied  by  Schmidt-Olt  et  al.  (1988). 
The  value  Dj-1.5  has  recently  been  reported  for  the  metal 
oxides  produced  by  spray  precipitation-oxidation  in  the  seeded 
premixed  flame  experiments  of  Matsoukas  and  Friedlander 
(1991).  In  this  classification  scheme  complete  restructuring  (but 
w’iihout  true  fusion,  which  would  ultimately  obliterate  evidence 
of  the  primary  spheres)  is,  of  course,  equivalent  to  the  familiar 
limiting  case  Dy  =  3,  corresponding  to  radially  constant  prop¬ 
erties  and  a  sharply  defined  outer  radius  larger  than 

the  gyration  radius  by  the  familiar  factor  (5/2)''^  as,  for 
example,  in  a  small  virgin  charcoal  or  polycrystalline  graphite 
particle  (Dixon-Lewis  et  al.,  1991).  [Some  authors  (for  ex¬ 
ample,  Ulrich,  1984)  distinguish  betw-ecn  physically  bound  ag¬ 
gregates  (easily  broken  apart)  and  partially  fusion-bonded 
agglomerates,  but  this  particular  distinction  will  not  concern 
us  here.]  The  goal  of  the  present  analysis  can  then  be  concisely 
stated  as  follows:  For  a  specified  aggregate  {N,  Dj)  composed 
of  a  large  number  N  of  primary  spheres  of  (nearly)  uniform 
radius  /?,,  on  the  surface  of  which  a  first-order  rate  process 
occurs  with  reaction  probability  a,  in  an  environment  char¬ 
acterized  by  (he  mean-free  path  /,,  what  is  the  ratio  of  the 
effective  (accessible)  surface  area  to  the  sum  of  the  primary 
sphere  surface  areas?  We  seek  a  rational,  yet  convenient  set 
of  quantitative  relationships  which  wdll  permit  such  calcula¬ 
tions  over  the  spectrum  of  carrier  fluid  mean  free-paths  from 
values  so  small  that  continuum  (bulk)  diffusion  occurs  even 
with  the  gas-filled  “pores,”  to  values  large  enough  to  not  only 
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lead  to  Knudsen  diffusion  within  the  pores  but  to  also  reduce 
the  diffusional  resistance  of  the  external  fluid  layer  surround¬ 
ing  the  aggregate.  The  procedures  we  ultimately  recorhmend 
are  efficient  enough  to  be  used  to  predict  the  accessible  area 
and  transport  behavior  of  populations  of  such  aggregates  but 
such  applications  (Rosner  and  Tandon,  1993a)  are  beyond  the 
scope  of  this  article.  Some  of  our  methods  and  results  would 
also  be  applicable  to  microporous  catalyst  supports  in  fixed- 
bed  reactors,  to  which  Sheintuch  and  Brandon  (1989)  have 
recently  applied  similar  pseudo-continuum  methods.  These  au¬ 
thors  treated  steady  and  nonsteady  Dj>  2.4  cases  in  the  absence 
of  external  diffusion  limitations  or  noncontinuum  (gas  kinetic) 
effects.  In  the  formulation  which  follows  we  explicitly  deal 
with  a  suspended  quasi-spherical  aggregate  comprised  of  a 
large  number  (N)  of  primary  spheres,  and  develop  a  conven¬ 
ient  rational  correlation  for  our  effectiveness  factor  results  in 
terms  of  the  prevailing  reaction  probability  or,  Knudsen  number 
based  on  primary  sphere  diameter,  Knm^-,  and  the  number  of 
primary  spheres  in  the  aggregate  over  a  broad  range  of  observed 
fractal  dimensions,  Dj. 

We  present  a  self-consistent,  tractable  pseudo-continuum 
mathematical  model  for  computing  the  accessible  area  of  a 
large  quasi-spherical  fractal  aggregate,  including  the  radial 
distributions  of  effective  pseudo-homogeneous  rate  constant 
and  reactant  molecule  effecti\e  diffusivity.  An  efficient  finite- 
analytic  method  is  implemented  for  ilie  required  numerical 
integrations,  and  a  convenient  correlation  procedure  is  devel¬ 
oped  based  on  the  introduction  of  an  effective  Thiele  modulus 
and  the  well-known  relationship  valid  for  the  limiting  case 
Z).  =  3.  Results  for  three  particular  £)^-values  (2.5,  2.18,  1.8) 
in  both  Knudsen  number  limits  are  discussed  as  w-ell  as  those 
which  account  for  arbitrary  Kn  and  external  transport  limi¬ 
tations.  Tests  of  the  recommended  effective  Thiele  modulus 
correlation  procedure  for  are  included.  Our  principal  as¬ 
sumptions  and  approximations  are  briefly  defended,  and  sev¬ 
eral  validity  criteria  and  tractable  generalizations  of  current 
interest  are  described.  Finally,  our  recommended  predictive 
procedure  and  its  application  to  coagulation-aged  aggregate 
populations  are  summarized. 


Mathematical  Model 

Basic  assumptions 

To  capture  the  essential  features  of  simultaneous  diffusion 
and  surface  reaction  in  a  large  aggregate  comprised  of  many 
(jV»  1)  spherical  primary  particles,  our  idealized  mathemat¬ 
ical  model  is  based  on  the  following  underlying  assumptions: 

(Al)  The  aggregate  behaves  like  a  motionless  quasi-spher¬ 
ical  “porous”  granular  solid  with  position-dependent  area  per 
unit  volume,  a"  and  effective  diffusivity  for  reactant  mole¬ 
cules, 

{A2)  The  radial  dependence  of  the  effective  first-order  rate 
constant  and  effective  diffusivity  for  reactant  molecules  are 
dictated  by  the  aggregate  (fractal)  exponent  Dj  defined  by 
dlnN/dlnr  (see  the  second  section  for  details). 

(A3)  A  first-order  chemical  (or  physical)  process  occurs  on 
the  external  surface  of  each  dense  primary  spherical  panicle 
comprising  the  aggregate  and  negligible  area  is  lost  at  the 
primary  particle  contact  points.  [The  second  part  of  this  as¬ 
sumption  must  ultimately  be  abandoned  w'ith  the  growth  of 


strong  fusion  bonds  between  primary  particles  (see  Koch  and 
Friedlander,  1990,  for  a  recent  summary  of  this  rate  process), 
or  loss  of  area  as  a  result  of  extensive  capillary  condensation 
or  growth  (CVD)  from  the  vapor.] 

(A4)  Reagent  A,  of  small  size  compared  to  the  primary 
particles,  is  steadily  supplied  to  the  intrinsic  surfaces  of  the 
aggregate  by  radial  diffusion  through  the  gas-filled  tortuous 
pores  defined  by  the  void  space  between  the  primary  particles. 

(A5)  The  aggregate  may  be  approximated  as  isothermal, 
that  is,  each  primary  particle  is  at  nearly  the  same  steady-state 
temperature  regardless  of  its  position  within  the  aggregate  (see, 
also,  the  section  discussing  primary  sphere  surface  roughness). 

(A6)  When  the  gas  mean-free-path  is  negligible  compared 
to  radial  concentration  diffusion  of  reagent  also  occurs 
throughout  the  space  external  to  the  aggregate, 

where  it  is  not  appreciably  impeded  by  the  presence  of  primary 
particles. 

We  postpone  to  the  fourth  section  a  defence  of  these  ideal¬ 
izations,  and  provide  criteria  to  define  their  approximate  limits 
of  validity.  Several  will  be  seen  in  the  fifth  section  to  be  readily 
generalized  to  embrace  more  complex  situations  that  will  un¬ 
doubtedly  arise  in  particular  applications.  Tentative  secondary 
assumptions  introduced  to  illustrate  our  procedures  (for  ex¬ 
ample,  the  nature  of  the  tortuosity-void  fraction  relation  for 
Knudsen-  or  bulk-diffusion)  will  be  introduced,  as  required, 
to  proceed. 


Aggregate  structure  for  Dj<3 

For  an  aggregate  whose  N{r)  relation  is  of  the  power-law 
form:  N  -  Pir/R^)^-',  where  (i  h  a  dimensionless  constant  dis¬ 
cussed  below,  w'e  have  the  basic  relation: 


dN=^Dj 


dr 


(1) 


describing  the  (expected)  number  of  primary  particles  within 
a  shell  of  thickness  dr  located  at  radius  r  (see  Figure  2).  [This 
level  of  description,  in  w'hich  Eq.  1,  in  effect,  defines  an  en¬ 
semble-averaged  radial  distribution  function  for  these  spatially 
nonuniform  random  structures,  will  prove  adequate  to  our 
present  needs.  However,  in  developing  theories  of  the  restruc¬ 
turing  kinetics  of  such  aggregates  (see  last  section)  additional 
microstruciural  information  is  needed,  and  w'e  introduce  and 
track  the  evolution  of  the  pdf  of  angles  between  triplets  of 
contacting  primary  particles  (Cohen  and  Rosner,  1994).]  This, 
in  turn,  implies  that  the  local  solid  fraction,  (pir)-,  is  given 
by: 


1^0  f  f  r\ 


(2) 


We  now  note  that  when  Z)/~3,  const,  so  that  the  di¬ 

mensionless  factor  0  (sometimes  called  the  “filling  factor”) 
may  be  identified  with  the  value  of  the  solid  fraction  in  the 
uniform  porosity  limit;  that  is,  /3  =  <^i,n,=  1  which,  for  a 

random  loose  packing  of  uniform  size  impenetrable  spheres, 
will  be  about  0.6.  Indeed,  we  visualize  that,  in  general,  Eq.  1 
applies  only  outside  of  a  uniform  porosity  “core”  of  radius 
where: 
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\  ^  ^  ^lim/ 


(3) 


This  equation  indicates  that  the  core  radius  is  necessarily  of 
order  /?,,  that  is,  quite  small  on  the  scale  of  the  overall  ag¬ 
gregate  size  defined  more  sharply  below. 

Whereas  most  aggregate  research  has  focused  attention  of 
the  radius  of  gyration,  we  instead  introduce  the  effective 
outer  radius,  defined  by: 


R 


max  ~ 


1/2 


(4) 


which  is  seen  to  reduce  to  the  true  outer  radius  of  the  spherical 
aggregate  in  the  Dj—  3  (uniform  porosity)  limit.  is  defined 
such  that  the  actual  rotational  moment  of  inertia  of  the  sphere 
is  the  same  as  if  all  of  its  mass  were  concentrated  at  its  gyration 
radius.  Rogak  and  Flagan  (1990),  who  were  primarily  concerned 
with  estimating  aggregate  particle^  r/rog  in  the  creeping  flow 
continuum  limit,  introduce  an  outer  radius  similar  to  that  given 
in  Eq.  4  but  omitting  the  factor  (3/2)'  ’^.  In  view  of  the  nature 
of  the  Z9y—3  limit  we  recommend  Eq.  4  in  such  future  studies. 
Mulholland  et  al.  (1988)  suggested  that  the  collision  radius  of 
free-molecule  clusters  (with  Dy=2)  is  about  2I®ro  greater  than 
Note  that  1(3/2)I(2  +  2)/21  r’=  |(3/2)  (2))'''^=  1 .73.  Thus, 
the  effective  collision  radius  appears  to  be  somewhat  less  than 
which  is  expected  on  physical  (partial  penetration) 
grounds.)  Beyond  /•=  treat  the  domain  as  primary  par- 

ticle-“free,’’  at  least  insofar  as  perimeter-averaged  reagent  ra¬ 
dial  diffusion  is  concerned — that  is,  we  treat  r  =  as  though 
it  were  a  “sharp”  boundary  of  the  large  N  aggregate  (Figure 

Summarizing,  if  we  define  the  normalized  radial  variable 
ci.r  idealized  model  of  a  large  quasi-spherical  ag¬ 
gregate  containing  N  primary  particles  is  as  follows: 


6  =  «iim*0-4  for  (0<{'<f.«  1)  (5a) 


3  \/3| 


for  (fr<  f<  1) 


(5b) 


for(l<f<oo)  (5c) 

where  c  is  (he  familiar  void  fraction  or  porosity,  1  -  .  and 

/3  =  0.6,  Whereas  higher-order  microstructural  informa¬ 

tion  (see,  for  example,  Torquato,  1991,  or  Tassopoulos  and 
Rosner,  1991,  1992)  may  be  needed  to  increase  (he  accuracy 
of  our  Knudsen  diffusion  coefficient  estimates  (below),  for 
our  present  purposes  the  above  level  of  information  will  prove 
to  be  sufficient — as  discussed  in  the  next  two  subsections. 


Bulk  diffusion  and  Knudsen  diffusion  lam:  associated 
tortuosities 

For  reagent  A  diffusion  within  the  aggregate  void  space  we 
must  distinguish  between  two  limiting  cases  depending  upon 
the  ratio  of  the  mean-free-path  /^,  to  the  local  mean  pore 
diameter.  When  this  ratio  (the  local  pore  Knudsen  number)  is 
very  small  we  write: 

(6) 


where  is  the  ordinary  Pick  diffusion  coefficient  for  spe¬ 
cies  A  migration  within  the  prevailing  gas  mixture  and  T,(f)  is 
the  continuum-limit  tortuosity  of  the  local  pore  structure  with 
void  fraction  e.  Similarly,  in  the  opposite  limit  (A'npf„,»l) 
we  write: 


(7) 


where  the  Knudsen  diffusion  coefficient  for  species  / 
migration  down  a  single  straight  pore  of  uniform  diameter,  is 
given  (for  diffuse  scattering)  by: 


D;c=(l/3)c,l(2/3)U/(l-0!(2/?.)l  (8) 


lOlherwise  should  contain  an  additional  factor  (2-/)//, 
where  /  describes  the  fraction  of  scattering  events  which  are 
“diffuse”  (cf.  “specular”).]  Here  c^  is  the  mean  thermal  speed 
of  gas  molecule  A  and  the  bracketed  factor  is  the  mean  pore 
diameter  for  a  random  packing  of  impenetrable  spheres  of 
radius  /?,  and  void  fraction  r.  While  recent  research  (for  ex¬ 
ample,  Tassopoulos  and  Rosner,  1992;  Melkole  and  Jensen, 
1992)  has  indicated  that  the  Knudsen  transition  is  probably 
somew'hat  more  complex,  for  arbitrary  Kn^^,,  we  here  adopt 
the  familiar  additive  resistance  relation: 


l^A.rU  ~  f ' 


(9) 


(see,  for  example,  Butt,  1980),  which  correctly  recovers  the 
above  mentioned  limititig  formulae,  Eqs.  6  and  7.  It  remains 
for  us  to  specify  the  indicated  continuum  and  Knudsen  tor¬ 
tuosities  for  the  random  packing  of  impenetrable  spheres  en¬ 
countered  in  such  self-similar  quasi-spherical  aggregates  as 
depicted  in  Figures  1  and  2. 

Recent  research  has  shown  that  T,(e)  is  !nicros(riic(ure-in- 
sensifive  and  often  adequately  described  by  a  power-law  in  e 
(Tassopoulos  and  Rosner,  1992).  Imposing  the  limits  t, -(1) 
=  1  and  Tr^O.4)-  =  1.48  (Huizenga  and  Smith,  1986),  we  there¬ 
fore  tentatively  adopt  (irrespective  of  D,): 

T,-(e,  (10) 


Knudsen  regime  tortuosities  are  known  to  be  more 

microstructure-sensitive,  and  vary  more  nearly  linearly  wdth  the 
local  void  fraction  €  (Tassopoulos  and  Rosner,  1992).  [Elias- 
Kohav  et  al.  (1991)  have  recently  calculated  Lhe  tortuosity- 
porosity  relation  for  a  variety  of  two-dimensional  fractal  porous 
media  and  contrasted  them  to  the  behavior  of  fully  random 
objects.  In  our  future  work  we  plan  to  probe  t  for  computer¬ 
generated  three-dimensional  aggregates  (cf.  the  analogous  work 
of  Tassopoulos  and  Rosner  (1991)  on  planar  granular  deposits.] 
Imposing  the  limits:  1 )-  =  1 ,  and  -(0.4}  =  1 .86  (Olague  el 

al.,  1988),  we  therefore  tentatively  adopt  (irrespective  of  Dj): 


T;,-(€}=r  1.86- 1.433(c~  0.40)  (11) 

While  the  methods  developed  below  will  readily  allow  the  in- 
troductfon  of  more  complex  tortuosity  law's  should  they  prove 
necessary,  all  specific  illustrations  contained  below-  (results  and 
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discussion  section)  will  be  based  on  these  particular  relations 
for  the  determination  of 


R  max 

1/2 

k 

First-order  irreversible  sink  law  for  /ferr 
Consider  a  first-order  rate  process  occurring  on  the  external 
surface  of  each  primary  particle.  For  a  collection  of  such  spheres 
the  available  surface  area  per  unit  volume,  a'\  can  be  written: 

a’"  :=l4>ir\/R,  (12) 

where  4>(r)  is  the  above  mentioned  local  solid  fraction.  There¬ 
fore.  the  pseudo-homogeneous  rate  constant  will  be  taken 
to  be: 


keU  —  ^ 

From  a  kinetic  theory  viewpoint  A''n,  the  heterogeneous  rate 
constant  (velocity)  can  be  reexpressed. 

A;n  =  (l/4)ac,  (14) 

where  a  is  the  dimensionless  reaction  probability  (a<  1)  and 
is  again  the  mean  thermal  speed  of  species  A.  [This  type  of 
kinetic  law  can  also  be  used  to  describe  physical  condensation 
(with  condensation  probability  a)  or  evaporation ,  with  evap¬ 
oration  sublimation  coefficient  cy,  to  be  discussed  later.]  W  hile 
a  is  not  really  an  elementary  rate  constant  (Rosner,  1972).  for 
our  present  purposes  we  w'ill  treat  q  as  a  specifiable  parameter. 


Pseudo-homogeneous  reaction-diffusion  model: 
spherical  symmetry 

Using  the  above  information  it  is  possible  to  derive  a  simple 
ODE  for  the  steady-state  concentration  profile  n^ir)-  estab¬ 
lished  within  the  aggregate,  viewing  the  latter  as  a  porous 
sphere  with  known  and  -profiles.  For  ex¬ 

ample,  a  species  A  mass  balance  for  a  spherical  shell  of  volume 
47rr^Ar  located  at  radius  r  (after  division  by  47rr'2kr  and  passing 
to  the  limit  Ar-O)  yields  the  familiar  linear  second-order  ODE: 


1  £ 

r^  dr 


dn^ 

~d7 


k^uir}n,ir^ 


(15) 


Introducing  the  normalized  variables. 


(16a, b) 


It  is  helpful  to  view  as  the  ratio  of  the  aggregate  radius  /?ma* 
to  the  characteristic  penetration  depth  defined  by  the  denom- 
inaior  of  Eq.  19.  Alternatively,  the  parameter  appearing 
on  the  righthand  side  of  Eq.  17  may  be  regarded  as  the  relevant 
Damkohler  group— that  is.  the  ratio  of  the  diffusion  lime 
to  the  reaction  lime:  (Ros¬ 

ner,  1986,  1990,  Chapter  8).  Thus,  once  the  coefficient  func¬ 
tions  and  ki^^  and  the  parameter  4>  are  specified,  it  is 

straightforward  to  find  c-(f>  satisfying  the  ODE  (Eq.  17)  and 
the  “split”  boundary  conditions: 

ci\)-  =  \  and  {dc/d^)^.o^O  (20a, b) 

(See  sections  on  finite-analytic  method  and  inclusion  of  ex¬ 
ternal  boundary  layer  limitations.)  The  first  of  these  BCs  fol¬ 
lows  from  the  definition  of  c  (Eq.  16a);  the  second  follows 
from  the  nonsingular  behavior  of  the  medium  as  r-0,  Ohat 
is,  within  the  constant  porosity  inner  core),  not  exclusively 
“symmetry”  (see  Rosner,  1986,  1990,  Section  6.4.4). 


Definition  and  calculation  of  the  "internal  effective¬ 
ness  factor 

Solution  to  the  boundary  value  problem  (BVP)  defined  in 
the  previous  subsection  allow's  calculation  of  the  single  quantity 
of  greatest  interest,  namely,  the  actual  reaction  rate  within  the 
agaregaie  compared  to  what  it  would  have  been  had  the  reagent 
concentration  remained  undepleted  everywhere  within  the  ag¬ 
gregate.  This  quantity,  hereafter  written  7;.,,,  can  also  be  re¬ 
garded  as  the  ratio  of  the  effective  (or  accessible)  surface  area 
to  the  true  total  primary  particle  surface  area  within  the  ag¬ 
gregate  (that  is,  the  sum  of  the  primary  sphere  areas,  S■a^). 
For  this  reason  this  quantity  is  traditionally  referred  to  as  the 
internal  effectiveness  factor, 

From  its  definition,  we  readily  find  that  7j,n,  can  be  calculated 
from  cin  via: 

■....4  (11  .  3  <3i) 


Thus,  once  the  solution  is  obtained  (see  ne.xt  subsection) 
the  internal  effectiveness  factor  (or  normalized  accessible  area) 
is  readily  found  from  Eq.  21.  It  is  interesting  to  note  that  since 
the  core  radius  makes  a  negligible  contribution  to  this  indicated 
integral,  Eqs.  2,  13,  18,  and  21  imply  that. 


then  Eq.  15  can  be  rewritten  in  dimensionless  form: 

1  — 

involving  the  coefficient  functions: 

k^k'.UrWk^liR.u.)-  dSa.b) 

and  the  parameter  <l>,  the  familiar  Thiele  modulus,  given  ex- 
pliciily  by: 


dc 


=  ¥kc 


(17) 


w'hich  is  a  straightforw-ard  generalization  of  the  familiar  Z>y=  3 
(constant  porosity)  result  (cf.,  for  example,  Rosner,  1986, 
1990).  However,  in  developing  a  compact  correlation  for  all 
of  our  computed  results,  Eq.  21  w-ill  prove  to  be  more  illu¬ 
minating  (see  subsection  after  next). 


Finite-analytic  method 

While  most  transport  (boundary  value)  problems  do  not  have 
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analytic  solutions,  some,  like  ODE  (Eq.  17)  and  its  BCs,  admit 
analytic  solutions  if  the  coefficient  functions  are  assumed  to 
be  constant  (here  and  'n  ^uch  ^'toalions  an 

attractive  scheme  to  solve  the  problem  is  to  simply  divide  the 
domain  into  sublayers,  each  with  piecewise-consian,  local  coef_ 
ficients.  The  BVP  over  the  whole  domain  can  then  be  solved 
by  applying  appropriate  continuity  conditions  at  the  bound¬ 
aries  of  each  sublayer.  Because  we  use  local  analytic  solutions 
applicable  within  each  sublayer  this  variant  of  the  f'mte-dif- 
ference  method  has  been  assigned  the  name:  "finite-analytic 
method  (Chen  and  Chen,  1984;  Rosner  et  al.,  1987;  and  Ros- 
ner,  1986,  1990,  p.  484).  For  any  fixed  number  of  sublayers 
we  expect’ the  error  involved  to  be  considerably  less  than  that 
involved  in  direct  discretization  of  the  governing  differential 
equation(s).  The  finite-analytic  approach  thus  facilitates  rapid 
convergence  to  the  exact  solution,  allowing  division  of  the 
domain  into  fewer  sublayers.  What  follows  is  a  brief  account 
of  its  implementation  for  the  problem  at  hand. 

For  constant  D  and  k,  the  exact  solution  to  Eq.  17  takes 

the  form: 


c=  \A  exp  +  Bexp  A  -  afH/f 


where 

With  the  boundary  conditions 
C-Ci,  the  constants  A  and  B  are: 


C  =  c,  and  f=h 


Figure  3a.  Discretization  nomenclature 

lation  of  a  finite-analytic  method  (FAM)  to 
solve  numerically  the  present  radial  diffu- 
slon/pseudohomogeneous  reaction  model 
for  the  accessible  area  of  a  large,  quasi- 
spherical  mulliparticle  aggregate  ot  nonun¬ 
iform  porosity. 


/4  =  lc,r,-exp-(oh)--Cih  _  . 

•exp  -(ofi)  VI^-XP  ~ 

fi=lc,  ^-exp•(o^>-f2h•exp-(-ofl)l/ 

[exp-(-2<jh)--exp-(-2ohH  (2-*) 


rfc/rfr=  [A  exp  ia^}  (of-  1)-Bexp  i  -  of)  (of -i-  l)]/f' 

For  a  porous  medium  with  fractal  exponent  D;<3,  the  coef¬ 
ficient  functions  D  and  k  are  not  constant  and  depend  on  f 
through  the  local  porosity  e-(f)  of  the  medium,  as  discussed 
earlier  Thus,  to  solve  the  problem  of  diffusion  and  reaction 
in  the  "fractal"  aggregate  we  assume  that  the  porous  medium 
consists  of  a  central  core  with  constant  porosity  enveloped  by 
successive  concentric  annular  shells.  The  porosity  (and  there¬ 
fore  l(  and  D)  is  assumed  to  be  constant  within  each  of  these 
shells  but  differs  from  shell  to  shell  (in  accord  with  e  evaluated 
at  each  shell  arithmetic  mean  radius),  as  given  by  Eq.  5.  The 
analytic  solution  as  given  by  Eq.  23  holds  within  each  of  these 
shells  for  the  local  values  of  l<  and  D.  As  illustrated  in  Figure 
3  the  outer  region  has  been  divided  into  n  zones  (shells)  an 
n’-v  1  nodes.  Defining  c,  and  f,  as  dimensionless  concentration 
and  distance  from  center  at  node  /,  respectively,  (cf.  Figure 
3a)  the  following  condition  should  be  satisfied  at  each  node  i: 


Substituting  from  Eqs.  23  and  24  we  get  a  linear  algeb^ak 

relationship  between  c,_ ,,  Cj,  and  c, 4  i  ^Vriling  Eq.  or  i 


r 

Figure  3b.  Typical  radial  profiles  of  normalized  pseudo- 
homogeneous  reaction  rate  coefficient  (k), 
reagent  diffusion  coefficient  (D),  and  corre¬ 
sponding  steady-state  normalized  composi¬ 
tion  profile  (c). 

Logarithmic  ordinaic;  cast  shown:  D,=  2.5.  10  .  a  =  10  , 

and  Xnj/g,  =  10. 

3  n  we  have  n  -  1  linear  algebraic  equations  in  n  +  1  vari¬ 
ables  (c,. /=  1,  ..  n-h  1).  Thetwoother  required  equations  are: 

.  .  =  1  (26) 
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and  species  A 


flux  continuity  condition  at  f-  Tcorc* 


^(dc\ 


=  {of-] 


Inside  the  core  the  concentration  profile  is  again  given  by  Eq. 
23  with  the  following  boundary  conditions: 


C  =  C|  @  f  =  fco>t 

(arc/rff)j-.„  =  0 

From  Eqs.  23  and  28  we  obtain; 

c  =  r,sinh-('J>corcr)-/(fs>"h-(4>cou3)  (29) 

where  4>,o,c  =  '(?>oK/(A-f(,«.o,c/^'=fr.«.ou)  •  cnlvprl 

This  system  of  n  +  1  linear  equations  can  be  rapidly  oh ed 
for  the  concentration  profile  (see,  for  example.  Figure  3b)  using 
Thomas-  algorithm  for  tridiagonal  matrices.  The  internal  e  - 
fectiveness  factor.  can  then  be  calculated  as  Eq.  21.  or. 
alternatively; 


forward  to  calculate  both  7)„,  and  via  simple  explicit 

algebraic  formulae. 

Just  as  our  finite-analytic  strategy  (previous  subsection)  tor 
the  solution  of  the  BVP  defined  earlier  was  based  on  the  fa¬ 
miliar  analytic  behavior  of  the  constant  property  (Dj-l)  spe¬ 
cial  case,  our  proposed  correlation  for  r,m  wiH  be  based  on  the 
well-known  properties  of  the  function: 

=(3/x)((ianhx)-'-Jr-'l  (31) 

which  (with  x  =  <i>)  describes  r,,„,  in  the  Dy=  3  limit,  and  which 
behaves  like  1  -  (2/5)x'-t- . . .  for  small  x  and  3/x  fo'' J^rge  x. 

Inspection  of  the  ODE  (Eq.  17)  suggests  that  when  t  is 
sufficiently  smalt  (to  allow  “deep  penetration"  of  the 
the  principal  effect  of  the  D;<3  terms  (which  lead  to 
Din<\  and  will  be  to  cause  the  solution  to  be 

similar  to  the  D/=3  case  but  with  an  effective  value  of 
“stretched”  by  the  factor: 

1<A'>/(D)]'''  (32) 

where  ( >  is  some  appropriate  average  value  of  the  indicated 
quantity  over  the  domain  0<f£  1.  We  postpone  for  the  mo¬ 
ment  the  question  of  which  average  is  most  successful  and  will 
empirically  consider  below  members  of  the  subclass. 


(>j,>,^(i+<7).5''k-(r>-m 


where  the  indicated  integrals  are  replaced  by  their  equivalent 

^'Tn^a  typical  Dy=  1.8  numerical  calculation  we  used  n=  10 
sublayers,  having  determined  that  at  even  this  level  ° 
division  our  results  are  quite  insensitive 
Sparc  II  station  run  time  was  less  than  5  s.  As  a  ^ 

on  the  computer  program,  in  the  specia  case  D,=  3  we  reco  e 
the  well  known  constant  property  result  for  t,.,„  (see  the  next 
subsection  and  Figure  7). 

Correlation  strategy  for  rO;<3) 

Despite  the  computational  efficiency  of  the  finite-analy_t.c 
numerical  method  described  in  the  previous  subsection  for 
accurately  calculating  r,,„,  for  any  particular  combination  of 
the  parameters  N,  Dj,  Kn,  . . .  and  choice  o  laws 
and  it  would  be  desirable  to  have  available  a  compaa. 

National  correlation  formula  which  could  be  used  to  rapid  y 
predict  n  with  acceptable  accuracy  (say  <  3^/o  rms  error)  ove 
fhe  emiie  range  of  aggregate/environmental  " 

engineering  interest.  This  is  especially  important  for  soK  ng 
prLlems  involving  the  presence  of  populations  of  aggregates 

for  which  N  covers  a  “spectrum”  of  values  and  we  agai 

to  calculate  the  total  accessible  surface  area  per  unit  volum 
represented  by  this  population  (see  second  to  last  subsection 
and  Rosner  and  Tandon,  1993a).  At  the  outset,  ^^^Id 
comment  that  while  this  calculation  requires  an  ^ 

^  which  is  the  product  of  and  a  corresponding  factor, 
“T'drsci^ing  the  effect  of  an  external  diffusion  layer,  we 
show  in  the  next  section  that,  once  t,,„,  is  known  it  is  stra.ght- 


where  ^  is  A  or  D  and  g>0. 

In  the  large  4>  limit  we  know  that  the  reagent  penetration 
will  be  shallow  so  that  A  =  1  and  D  =  1  within  this  outer  “bound¬ 
ary  layer.”  Thus,  in  this  limit,  Eq.  21  reveals  that  the  solution 
to  our  BVP  must  be: 


=  3  =31<AO:4>r 

,p»\  C  J 

in  which  the  relevant  effective-<t.  is  seen  to  be 
the  factor  <A->,.  Summarizing,  we  expect  that  r)3-(4>  „k  will  be 
an  acceptable  approximation  to  i,,„,  provided  <t>,„  is  chosen  such 
that  the  stretching  factor  l<A->/<D>l''‘ '.-rFor 

4.  and  the  stretching  factor  <A>j  is  applied  for  large  4>.  For 

these  reasons  we  simply  investigate  the  choices. 


4>,.,i-(‘l')-  ® 


For  each  such  choice  we  compared  values  of  n,,,,  (correlation) 
with  T,,„,  (FAM).  We  found  that  all  of  our  D;<3 
were  quite  acceptable  for  the  choice  9  =  1  •  This  is  readily  dem¬ 
onstrated  in  the  “correlation”  plots  Figures  4,  5  and  6  for 
cases  involving  continuum  (bulk)  pore  diffusion  (F‘f 
Knudsen  pore  diffusion  (Figure  5).  and  transitional  knudsen 

numbers  (Figure  6). 

One  sees  that  the  rms  errors  are  only  about  3  .0.  ’^e 

acceptably  small  local  errors  peaking  (ca.  4«?o)  at 
n  -values  (near  0.8).  Thus,  to  rapidly  predict  r),„,  for  larg 
aggregates  to  within  about  ca.  3%  error  one  merely  uses  the 
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Figure  4.  Test  of  the  effective  Thiele  modulus  correla¬ 
tion  (Eqs.  31  and  35)  for  the  internal  effective¬ 
ness  factor  nin,  of  large  (W>10*)  fractal 
aggregates  in  the  continuum  (bulk-)  pore  dif¬ 
fusion  limit. 

For  unshaded  points  10;  tot  shaded  points:  I. 


well-known  tjj-tx)-  function  with  a  suitably  calculated  (nonliii- 
early  stretched)  effective  Thiele  modulus,  ♦at  (Eq.  35)  with 
7=1.  [We  suggest  that  this  rational  correlation  strategy  can 
be  implemented  for  a  wide  class  of  engineering  problems.  Ear- 


Figure  5.  Test  of  the  effective  Thiele  modulus  correla¬ 
tion  (Eqs.  31  and  35)  for  the  internal  effective¬ 
ness  factor  ijini  of  large  (A/>10’)  fractal 
aggregates  in  the  Knudsen  pore  diffusion 
limit. 


Figure  6.  Test  of  the  effective  Thiele  modulus  correla¬ 
tion  (Eqs.  31  and  35)  for  the  internal  effective¬ 
ness  factor  Tim,  of  large  (WalO^)  fractal 
aggregates  in  the  transition  pore  diffusion  re¬ 
gime;  Knjn,  =  10. 


Her,  closely  related  examples  from  the  author’s  research  (DER) 
may  be  found  in  Israel  and  Rosjier  (1983)  and  Rosner  and 
Feng  (1974).]  Since,  in  this  case  (k)„  can  be  evaluated  in  closed 
form  (for  any  q)  this  recommended  procedure  only  requires 
the  calculation  of  <D>,,  which  is  a  trivial  numerical  integration 
for  the  present  class  of  tortuosity  laws  (second  section),  or  any 
other  choice  of  tortuosity  laws  deemed  (via  future  research  on 
aggregate  microstruclure)  to  be  more  appropriate  (fourth  sec¬ 
tion). 

In  closing  this  section,  we  remark  that  the  acceptably  small 
local  errors  of  the  present  recommended  correlation  appear  to 
be  systematically  positive  (see  Figures  4,  5,  6).  This  implies 
that  a  slight  modification  of  the  exponents  appearing  in  Eq. 
35  could  easily  reduce  them  further.  However,  the  present 
simple  version  (the  first  we  investigated!)  is  acceptable  for  our 
present  purposes.  We  remark  that  this  procedure  should  also 
be  useful  for  estimating  ii.n,  for  ordinary  catalyst  pellets  (Df  =  3) 
which,  because  of  their  manufacture  (pelleting,  impregnation, 
and  so  on),  hz'/e  spatially  nonuniform  D,a..m  srid  A  ,,,  properties. 

Results  and  Discussion 

Inclusion  of  external  boundary  layer  limitations 

If  the  gas  mean-free-path  is  small  on  the  scale  of  (Eq. 
4),  then,  even  though  the  aggregate  may  be  iVi/er/io//j' accessible 
(as  when  D/=  1.8  and  Knudsen  diffusion  prevails  within  the 
aggregate)  the  steady-state  reagent  concentration  established 
at  r  =  will  be  systematically  less  than  that  in  the  ambient, 

contributing  to  an  apparent  reduction  in  accessible  surface 
area.  This  effect  of  the  “resistance  in  series”  represented  by 
the  external  gas  boundary  layer  is  easily  included  in  our  analysis 
via  introduction  of  the  additional  factor; 
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.00 


(36) 

v.hich  is  explicitly  calculated  below.  The  overall  effect  of  in¬ 
ternal  and  external  diffusional  resistance  on  the  normalized 
accessible  surface  area  of  the  aggregate  will  then  be  given  by 
the  product: 


where  Kn.n,  is  the  Knudsen  number  based  on  primary  sphere 
diameter ’2/?,,  a  quantity  specified  explicitly  in  the  aggregate 
transport  calculations  illustrated  below'.  Thus: 

A'/7  = J 

-  ) "  ‘  =  •  (AV/3) "  (43) 


^o^^■rall  (37) 

where  the  calculation  and  correlation  of  have  already  been 
dealt  with  earlier. 

Our  idealized  mathematical  model  to  account  for  the  ex¬ 
ternal  diffusion  layer  makes  use  of  earlier  results  for  the  di¬ 
mensionless  (Npsselt-Sherwood)  mass-transfer  coefficient  for 
diffusional  transfer  to/from  a  motionless,  isolated  solid  sphere, 
viz. 


iVu„,  =  2- 1 1 +/r/-{AV/>r'  (38) 


where  Kn-  1/{2R^^^).  Note  that  /c/(0)-0  in  order  to  recover 
the  familiar  continuum  limit  value  2,  but  fctico)-^  oo  to  elim¬ 
inate  this  diffusional  resistance  in  the  free-molecule  limit.  Sev¬ 
eral  authors  have  calculated  and/or  measured  the  important 
function /r/^A'/jf  appearing  in  Eq.  38.  For  our  present  pur¬ 
poses  we  adopt  the  following  corrected  form  of  the  Fuchs- 
Sutugin  approximation  (Friedlander,  1977); 


Jct{Kn)  =  IKn  ■ 


.33 1.42(2AV?)' 
\  ^{2Kn)  - ' 


(39) 


(The  coefficient  of  (2A'/t)‘'  in  the  numerator  of  Eq.  39  has 
been  increased  (from  0.71)  to  force  agreement  with  independ¬ 
ent  Kn«  1  results.]  Then  a  steady-state  analysis  in  which  the 
reagent  A  flow  rate  into  the  quasi-spherical  aggregate  is  bal¬ 
anced  against  the  steady-state  rate  of  its  consumption  within 
the  aggregate  yields  the  explicit  algebraic  result: 


-  /  1  + 


^  A  •  77, r, I 


(40) 


where  (  =47r/?f)  is  the  surface  area  of  a  primary  spherule, 
For  the  illustrative  calculations  given  below'  we  will  assume 
that  reagent  A  has  a  molecular  mass  and  collision  cross-section 
similar  10  that  of  the  background  gas  mixture,  so  that  4  =  /^. 
We  w'ill  also  make  use  of  the  corresponding  molecular  theory 
estimate: 


(41) 


These  (unnecessary)  assumptions  will  simply  reduce  the  num¬ 
ber  of  parameters  that  must  be  specified  in  the  examples  below', 
and  thereby  help  us  focus  on  the  essential  trends.  With  these 
simplifications  we  find  that  can  be  expressed; 


.^int 


Nu„,  •  Kn^ 


- 1 


(42) 


which  fixes  Kn  (cf.  Eq.  39)  for  every  choice  of  aggregate  size 
N  and  Knudsen  number  based  on  primary  sphere  diameter. 

Illustrative  results  and  implications 

Before  examining  representative  results  for  it  is  inter¬ 
esting  to  examine  the  sensitivity  of  77,,,,  to  the  structural  exponent 
Dj  and  aggregate  size  N.  One  such  set  of  results  is  shown  in 
Figure  7,  which  reveals  that  even  for  a  reaction  probability  of 
unity,  and  for  aggregates  containing  very  many  primary  par¬ 
ticles,  if  Df  is  less  than  2  and  Kn^^^  is  large  enough  (here 
-  10)  essentially  the  entire  aggregate  area  is  “accessible” 
-  !)•  [Indeed,  for  Dj<2  one  finds  that  for  Knudsen  pore 
diffusion  \\\t  penetration  depth  Jl’"  in¬ 

creases  more  rapidly  with  A/ than  itself,  that  is,  the  relevant 
Thiele  modulus  ^  actually  decreases  w'ith  increasing  AM  How¬ 
ever,  even  for  Dj<2  this  remarkable  behavior  does  not  occur 
in  the  continuum  pore  diffusion  limit  (A^/7p^,,<  1).]  However, 
for  a  given  A'-value  (say,  N=  10^)  one  sees  that  77,,,,  falls  off 
rapidly  at  higher  values  of  the  exponent  Dj.  Thus,  if  for  any 
reason  restructuring  of  an  initially  Dj=  1.8  aggregate  containing 
A' primary  particles  occurred,  77, w’ould  drop  accordingly,  even 
if  primary  sphere  fusion  did  not  occur. 

In  Figure  8,  which  serves  as  one  check  on  our  general  com¬ 
puter  program  (incorporating  the  pore  diffusion  transition  be- 
tw'een  continuum  and  Knudsen),  we  show  how  77,,,,  relates  to 
its  (separately  computed)  continuum-  and  Knudsen-pore  dif¬ 
fusion  asymptotes.  Note  that,  in  the  “small”  N  domain,  be¬ 
cause  of  the  reduced  average  pore  size,  one  approaches  the 
Knudsen  diffusion  asymptote,  whereas  for  N»  1  more  of  the 
aggregate  pores  fall  in  the  continuum  (bulk)  diffusion  regime. 


NUMBER  OF  PRIMARY  SPHERES.  N 

Figure  7.  Predicted  dependence  of  the  Internal  effec¬ 
tiveness  factor,  on  a  aggregate  size  N  and 
fractal  dimension  Df  for  a=  10“’;  Kn2R^  =  10. 
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LOG  N 

10 

Figure  8.  Nature  of  the  transition  from  obstructed 
continuum  (bulk  )  diffusion  to  Knudsen  pore 
diffusion  for  large  quasi  spherlcal  aggregates 
(D,=  2.5,  Kn2p,  =  10,  a=1). 


This  particular  graph  applies  to  the  parameter  combination: 
Df=2.5y  0=1,  Kn2K,=  \0.  However,  as  commented  above, 
even  for  Dy=  K8  low  values  oT  77, result  at  small  enough 
ambient  mean  free  paths  (high  enough  pressures). 

For  the  same  value  of  Kn^R^  Figures  9  and  10  reveal  the 
sensitivity  of  the  overall  normalized  accessible  surface  area, 
T?ovpiaii»  to  aggregate  size  N  and  reaction  (or  “trapping”)  prob« 
ability  a,  for  the  important  special  cases  Df~  1.8  (open  ag¬ 
gregates)  and  Dj-  2.5  (more  compact  aggregates).  Wc  remark 
that  the  exponent  Dj~  1 .8  ( =t  0. 1 )  has  been  observed  or  inferred 
for  aggregate  smokes  in  many  combustion  environments  (see 


NUMBER  OF  PRIMARY  SPHERES,  N 

Figure  10.  Predicted  dependence  of  the  overall  effec* 
tiveness  factor,  on  aggregate  size  N 

and  reaction  probability  «  for  D,=  2.5 

the  useful  summaries  of  Megaridis  and  Dobbins  (1990)  and 
Puri  et  al.  (1993).  [Megaridis  and  Dobbins  (1990)  have  also 
called  attention  to  the  narrow  size  range  (spread)  of  the  primary 
particle  diameters  in  particular  laminar  flow  environments. 
Primary  particle  diameters  reported  are  often  of  the  order  of 
tens  of  nanometers,  which  accounts  for  the  very  high  specific 
surface  areas  of  these  materials.)  The  value  D^  =  2.5  was  re¬ 
cently  observed  for  aggregated  metal  oxide  smokes  produced 
by  spray  pyrolysis-oxidation  of  aqueous  metal  salts  in  laminar, 
atmospheric  pressure  premixed  flames  (Matsoukas  and  Fried- 
lander,  1991).  Monte  Carlo  simulation  studies  in  3  dimension.s 
(Meakin,  1983,  1984)  reveal  that  the  Dy-valucs  near  1.8  result 
from  cluster-cluster  agglomeration  w-hcreas  Devalues  near  2.5 


NUMBER  OF  PRIMARY  SPHERES.  N 

Figure  9.  Predicted  dependence  of  the  overall  effec¬ 
tiveness  factor,  on  aggregate  size  N  and 
reaction  probability  a  for  Df^  1.8  (Kn2fl,  =  10). 


NUMBER  OF  PRIMARY  SPHERES,  N 

Figure  11.  Predicted  dependence  of  the  overall  effec¬ 
tiveness  factor,  on  aggregate  size  N 
and  gas  mean  free  path  (In  multiples  of  2R^) 
for  a  reaction  probability  a  of  10~^  1.8). 
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result  from  cluster-monomer  aggregation  in  the  absence  of 
further  restructuring  (Figure  1).  Whatever  the  prevailing  Dr 
value,  our  procedure  is  seen  to  facilitate  rational  predictions 
of  the  accessible  surface  area  of  large  quasi-spherical  aggre¬ 
gates  as  a  function  of  both  reaction  probability  and  Knudsen 
number. 

Whereas  some  authors  have  stated  or  implied  that  aggregates 
with  Dj<2  are  alw'ays  “totally  accessible”  this  will  clearly  not 
be  the  case  at  pressure  levels  large  enough  to  cause  the  mean- 
free-path  to  be  comparable  with  the  average  pore  diameter 
|(2/3)[t/(l -f)](27?,)  1 ,  or,  certainly,  the  primary  particle  di¬ 
ameter.  Moreover,  unless  /,  is  large  compared  to  the 

external  boundary  layer  limitation  is  not  escapable.  This  is 
clearly  illustrated  in  Figure  II  for  the  case  of  Dj-  1.8  with  a 
reaction  probability  of  10'*.  We  remark  that  reaction  prob¬ 
abilities  of  approximately  this  order  of  magnitude  have  been 
reported  for  the  oxidation  of  suspended  carbonaceous  soot  by 
OH(g)  or  0(g)  (Neoh  ei  al.,  1981;  Roth  et  al.,  1990).  In  con¬ 
trast,  inferred  reaction  probabilities  for  the  growth  of  “young” 
soot  via  C2H2(g)  impacts  are  of  the  order  of  10'^  (Harris  and 
Weiner,  1984).  Our  approach  and  results  may  therefore  be  of 
special  interest  in  high-pressure  devices--such  as  modern  air¬ 
craft  gas  turbine  combustors  (which  operate  at  ca.  30  atm, 
Correa  and  Shyy,  1987),  or  flame-produced  Ti02(s)  (chloride- 
process)  pigment  reactors  (Ulrich,  1984).  In  this  respect,  re¬ 
examination  of  Eq.*41  used  above  indicates  that,  provided 
assumptions  1-6  (second  section)  are  satisfied,  our  present 
results  could  also  be  used  to  estimate  the  accessible  area  of 
such  aggregates  in  supercritical  vapor.s  (Mohamed  et  al.,  1989) 
or  true  liquids  (Matijevic,  1981;  Zukoski  et  al.,  1990}  by  setting 
A^w,„  =  2,  and  reinterpreting  as  a  dimensionless  solute 

diffusion  coefficient  (inverse  Peclei  number): 

J„p/ 1  (( 1  /3)r, )  •  ( 2/? , )  |  (44) 

where  ( A4  mix)t»p  js  the  experimental  solute  A  Pick  (or  Browm- 
ian-)  diffusion  coefficient  in  the  pre\ailing  fluid,  and  the  mean 
molecular  speed  is  formally  calculated  using  gas  kinetic 
theory.  In  this  respect  the  effective  reaction  probability, 
should  be  calculated  from: 

aeff=(^'^)cxp/l(l/4)pj  (45) 

Clearly,  how'ever,  these  results  would  have  to  be  generalized 
to  deal  with  (“hindered”  diffusion)  situations  in  which  the 
effective  size  of  a  (macro-)  solute  A  is  not  negligible  compared 
ioiheprevaiiingspherulediameter(cf.  A4).  indeed,  Hagenlock 
and  Friedlander  (1989)  have  carried  out  Monte  Carlo  simu¬ 
lations  of  finite  size  gas  molecules  striking  Dj^2,5  aggregates 
of  spheres  with  N<  i0\  keeping  track  of  collision  frequency 
(assuming  diffuse  reflection)  as  a  function  of  aggregate  size 
and  gaseous  Knudsen  number  /^/(2y?,).  Remarkably,  they  re¬ 
port  steady-state  collision  frequencies  less  than  that  expected 
for  isolated  spheres  even  in  the  nonreactive  limit  (a~0)  and 
point  molecule  limit  (a^/(2/?, )  —  0),  a  paradoxical  result  evi¬ 
dently  at  variance  with  the  predictions  of  our  present  pseudo¬ 
continuum  model  (for  which  1  when  a  — Oand  a^/(2/?,)— 0 
for  any  fixed  N,  Kn). 

Discussion 

Equations  15-43  w'ill  clearly  permit  rapid  estimates  for  many 


parameter  combinations  of  particular  interest  to  the  reader. 
Moreover,  as  mentioned  in  the  fifth  section  and  will  be  ex¬ 
plicitly  illustrated  in  Rosner  and  Tandon  (1993a),  they  are 
readily  applied  to  predict  the  accessible  area  of,  say,  log-normal 
populations  oi  aggregates  frequently  observed  in  coagulation- 
aged  systems  (Matsoukas  and  Friedlander,  1991). 

It  is  noteworthy  that,  whereas  Meakin  and  Whitten  (1983), 
using  Monte  Carlo  simulation  methods,  have  implied  that  the 
accessible  area  (or  interface)  of  an  aggregate  will  scale  as  a 
simple  power  (near  0.74  in  3  dimensions)  of  the  total  aggregate 
size,  our  pseudo-continuum  relations,  which  also  include  ex¬ 
ternal  boundary  layer  effects,  reveal  that  the  actual  relationship 
will  not  be  so  simple.  This  is  readily  seen  by  examining  the 
range  of  values  of  1  -i-  (d in N)  revealed  in  Figures 
9,  10,  and  1 1.  (Note  that  the  logarithmic  coordinates  of  these 
figures  allow  direct  visualizaiion/computation  of  the  required 
local  logarithmic  derivative.) 

It  is  also  interesting  to  note  that  the  present  approach  can 
shed  useful  light  on  the  meaning  of  recent  experimental  meas¬ 
urements  suggested  to  characterize  aggregates  in  situ.  One  such 
measurement  recommended  and  illustrated  by  Schmidt-Oii 
(1988)  is  the  so-called  attachment  coefficient  for  aggregate 
scavenging  of,  say,  radioactive  ''‘Pb(g)-atoms.  In  our  notation 
this  quantity  (called  /3  by  Schmidi-Ott)  is  equivalent  to: 

O  ~  (  1  /4)aC^  a^'  N  ^overatl  (46) 

where  the  important  product  may  be  regarded  as  the 

accessible  (or  “exposed”)  number  of  primary  particles,  and, 
in  this  case,  >4  =^'‘P/}(g).  In  the  near  free-molecular  limit 
(A';7»  1)  this  coefficient  has  been  shown  to  be  approximately 
proportionalio  photoelectron yield (QunschtT  2Ln(\  Schmidi- 
Ott,  I985)of  small  (N< 40,  £)/  — 2.18, /?|  =  7.5  nm)  aggregates 
of  silver  illuminated  by  uv  radiation  (X  =  0-(200  nm}.  [This 
is  already  an  example  for  w'hich  should  be  distinguished 
from  the  gas  mixture  mean-free-path,  and  Eq.  41  for 
suitably  generalized  (Rosner,  1986,  1990).]  However,  inspec¬ 
tion  of  Eq.  42  and  the  laws  of  electron  scattering  in  gases 
suggests  that  this  linearity  is  not  likely  to  be  preserved  in  the 
continuum  (high-pressure)  limit. 

It  is  also  noteworthy  that  there  may  be  a  sort  of  “Reynolds 
analogy”  relation  (Rosner,  1986)  between  the  mass-transfer 
(scavenging)  behavior  of  aggregates  and  their  momentum 
transfer  behavior.  Indeed,  a  near  equality  between  the  effective 
scavenging  radius  and  the  mobility  equivalent  radius  has  been 
reported  and  discussed  by  Meakin  et  al.  (1989),  Schmidi-Oti 
et  al.  (1990),  and  Rogak  et  al.  (1991)  (who  report  agreement 
within  ca.  15%  for  Kn  <  0(1)).  This  implies  that  our  scavenging 
rate  formalism  for  aggregates,  at  least  in  the  high  a  limit, 
provides  some  insight  as  to  the  mobility  and  hence  Brownian 
diffusiviiy,  of  such  aggregates  over  a  wide  range  of  sizes  and 
A'/72k, -values.  In  view  of  the  simplicity  and  computational  ef¬ 
ficiency  of  our  methods  these  implications  certainly  warrant 
future  investigation. 

Defence  of  Approximations,  Generalizations 

It  would  be  prudent  to  briefly  reconsider  some  of  our  fun¬ 
damental  approximations,  indicating  their  expected  domain  of 
validity  and  suggesting  possible  generalizations  of  likely  future 
interest. 
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Continuum  approximation  for  small  N  and  Dj 

Even  though  the  simple  scaling  law:  d In  N/d /at  r=  const. ^ 
has  been  found  to  apply  down  to  remarkably  small  aggregate 
sizes  (Schmidt-Olt  et  al.,  1988,  1990;  Megaridis  and  Dobbins, 
1990)  clearly,  our  pseudo-continuum  mathematical  model  (in 
which  each  aggregate  is  treated  as  a  quasi-sphcrical  porous 
object)  should  not  be  expected  to  apply  to  aggregates  con¬ 
taining,  say,  less  than  about  30  primary  particles.  While  the 
extent  to  which  such  pseudo-continuum  results  match  small  N 
detailed  particle-level  diffusional  simulations  (Rosner  cl  al., 
1991;  Mackowski,  1994)  remains  to  be  investigated,  we  should 
therefore  not  expect  our  present  results  to  be  quantitatively 
accurate  for  aggregate  sizes  much  below  about  I0^  especially 
when  Dj  is  less  than  2.  jThe  situation  here  is  quite  similar  to 
our  earlier  theoretical  studies  of  the  evaporative  combustion 
of  a  “cloud”  containing  N  dispersed  fuel  droplets  (Labowsky 
and  Rosner,  1978).  In  anticipation  of  relevant  “small  A/”  re¬ 
sults  we  have  deliberately  started  our  graphs  al  10,  in  order 
to  embrace  some  of  the  expected  “matching”  domain.] 

Tortuosity-porosity  laws  over  the  spectrum  of  pore 
Knudsen  numbers 

Current  research  (extensions  of  Tassopoulos  and  Rosner, 
1992;  Melkote  and  Jensen,  1992;  Elias-Kohnv  ct  al.,  1991)  may 
reveal  that. our  presently  used  tortuosity  laws  are  not  suffi¬ 
ciently  accurate,  especi  illy  in  the  Knudsen  pore  diffusion  limit 
(w'here  higher-order  microstructural  information  characteriz¬ 
ing  the  likely  primary  sphere  positions  may  need  to  be  added) 
or  in  the  transition  regime  (where  the  pore  Knudsen  number 
dependence  of  the  tortuosity  may  not  quite  be  monolonic  and, 
accordingly,  depart  loo  much  from  the  behavior  tentatively 
assumed  here).  W'hatever  the  outcome  of  such  studies,  how¬ 
ever,  the  present  pseudo-continuum  procedures  should  remain 
useful  to  (re-)compulc  and  correlate  the  accessible  surface  area 
of  such  aggregates  as  a  fraction  of  the  total  area  of  the  primary 
spheres  present.  Indeed,  in  the  near  future  it  would  be  inter¬ 
esting  to  compare  the  present  estimates  with  much  more 
computationally-intcnsive  Brownian  simulations  on  computer- 
generated  fractal  aggregates  in  3  dimensions. 

In  closing  this  section  it  should  be  remarked  that  unusual 
Knudsen  diffusion  conditions  would  occur  deep  inside  such 
an  aggregate  since,  as  one  approaches  the  “core”  one  can 
show  that  the  local  characteristic  radial  length  over  which  the 
porosity  changes  by  an  appreciable  fractional  amount,  that  is, 
{dlnt/dr)'^  will  no  longer  be  large  on  the  scale  of  the  local 
average  pore  diameter!  We  should  also  comment  that  the  limit 
a— 1,  w'hile  formally  included  here,  should  also  be  expected 
to  be  singular  for  Knudsen  diffusion  (when  of=  I  no  reagent 
molecule  survives  one  collision  with  the  pore  w^all!).  This  im¬ 
plies  that  corrections  to  Eq.  8  should  be  introduced  in  the  high 
reaction  probability  limit  (Verhoff  and  Sireider,  1971). 

Nonspherical  aggregate  shapes 

It  has  been  reported  that  “small”  aggregate-aggregate  en¬ 
counters  lead  to  larger  aggregates  for  which  the  gyration  radii 
are  not  the  same  in  all  directions.  Bold  and  Jullien  (1986)  and 
Hentschel  (1984)  report  an  asymptotic  gyration  radius  ratio  of 
ca.  2.  Moreover,  electron  micrographs  of  soots  sampled  from 
diverse  systems  usually  reveal  a  wide  variety  of  nonspherical 


aggregate  morphologies,  often  somewhat  stringy  in  character 
(Ulrich,  1984;  Lahaye  and  Prado,  1981).  While  there  is  ample 
precedent  for  treating  the  effectiveness  of  nonspherical  porous 
catalyst  pellets  (with  Dy=3)  using  the  notion  of  an  equivalent 
sphere  of  effective  radius: 

~  3  ■  (  / A  )pulcr  rnvciopt  (47) 

(See  Aris,  1957).  If  increased  accuracy  were  needed  the  present 
model  could  readily  be  generalized  to  deal  with,  say,  prolate 
spheroidal  porous  media  which  are  fractal-like  in  each  principal 
direction.  Moreover,  external  diffusion  layer  transport  could 
be  dealt  with  using  the  dimensionless  transfer  coefficient  for 
such  spheroids  (Yovanovich.  1987)  but  preferably  based  on 
the  characteristic  length  These  generalizations 

would  allow  for  the  fact  that  the  characteristic  Brownian  ro- 
tation  time  of  such  large  aggregates,  for  example,  (6Dx,p,) 
would  be  much  larger  than  the  characteristic  time  (  {A.Jtt)/ 
for  reactant  molecule  diffusion  across  the  exter¬ 
nal  boundary  layer,  especially  at  small  Knudsen  numbers  based 
on  the  length  (/1„,/7r)''^  Thus,  the  often-quoted  notion  that 
Browmian  rotation  would  ensure  the  reasonableness  of  a  quasi- 
spherical  aggregate  approximation  even  for  large  elongated 
aggregates  is  probably  illusory.  How-ever,  we  leave  such  gen¬ 
eralizations  to  the  future,  noting  only  that  some  of  the  features 
of  porous  suspended  particles  which  are  distributed  with  re¬ 
spect  to  both  size  N  and  aspect  ratio  /R  are  treated  in  Rosner 
and  Tandon  (1994). 

Kinetic  layvs  other  than  first-order  irreversible 

As  mentioned  earlier,  the  present  treatment  for  and  T^o'cian 
can  be  generalized  to  apply  equally  well  to  simple  physical 
phenomena  such  as  reversible  condensation  or  evaporation  in 
w'hich  the  interfacial  rate  process  (sink  strength)  is  proportional 
to, a  \n^-n^,^iT,.\  p}]  where  a  is  now  the  condensation  or 
evaporation  coefficient.  How'ever,  in  some  systems  the  dt- 
pendence  of  or  and  \  p}  on  the  radius  (of  curvature) 

of  the  primary  particles  may  have  to  be  explicitly  taken  into 
account.  More  complex  nonlinear  kinetic  laws  w'ill  require 
separate  treatment  and  will  inevitably  reveal  a  dependence  on 
additional  dimensionless  parameters,  as  is  already  familiar  in 
the  field  of  heterogeneous  catalysis  (Aris,  1975;  Satterfield. 
1970;  Froment  and  Bischoff,  1979;  Engasscr  and  Horvath, 
1973). 

Primary  sphere  surface  roughness  and/or  internal 
porosity 

Using  X-ray  scattering  or  adsorption  techniques  some  ag¬ 
gregated  soot  systems  have  been  reported  to  exhibit  significant 
primary  particle  surface  roughness  and/or  porosity  within  the 
primary  spherules  (Lahaye  and  Prado,  1981 ;  Hurd  ct  al.,  1987). 
If  necessary,  such  behavior  could  be  modeled  by  imagining 
the  primary  particles  to  consist  of  an  outer  microporous  layer 
covering  an  impermeable  core,  or  containing  micropores 
throughout.  Indeed,  in  the  spirit  of  self-similarity  (Feder ,  1988; 
Dimotakis,  1991),  what  we  have  called  the  primary  particles 
could  themselves  be  viewed  as  fractal  microspheres,  thereby 
making  direct  use  (on  a  smaller  scale)  of  the  concepts  described 
here.  From  this  vantage  point  our  present  treatment  (As- 
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sumption  A3)  is  therefore  seen  to  be  equivalent  to  assuming 
that  the  primary  spheres  are  themselves  characterized  by  an 
effective  primary  sphere  Thiele  modulus  so  large  that  their 
internal  effectiveness  factors  are  very  small,  corresponding  to 
reaction  or  trapping  on  their  outer  surfaces. 

Nonisothermal  behavior 

Since  most  of  the  interfacial  processes  mentioned  above 
(physical  condensalion/evaporation,  reactive  growth  from  the 
vapor  phase,  oxidation  (gasification)  by  one  or  more  reagents 
in  the  vapor  phase)  are  not  ergoneuiral;  in  general  there  will 
be  a  tendency  for  the  aggregate  to  overheat  or  undercoo!  rel¬ 
ative  to  the  surrounding  gas  mixture.  Moreover,  if  the  effective 
Biot  number,  Bif,,  for  radial  energy  transfer  defined  by: 

5/XV^a„,trf)‘(>V/v,/2)  (48) 

(Rosner,  1986,  1990)  is  not  sufficiently  small,  then  temperature 
nonuniformities  within  the  aggregate  may  play  a  role  in  de¬ 
termining  both  17, nt  and  r/ovcraii.  especially  when  the  reaction 
probability  a  is  low  and  sharply  (Arrhenius-)  temperature  de¬ 
pendent,  and  the  mole  fraction  of  the  reagent  is  not  very  small. 
(In  this  case  it  might  also  be  necessary  to  account  for  the 
tendency,  due  to  radiative  energy  loss,  for  the  larger  aggregates 
to  be  systematically  cooler  than  the  smaller  ones  (Mackowski 
et  al.,  1991,  1994;  Rosner  et  al.,  1992).] 

At  the  moment,  not  enough  is  known  about  the  local  ef¬ 
fective  radial  thermal  conductivity,  within  such  aggre¬ 

gates,  including  the  Rayleigh-limit  radiative  contribution,  to 
realistically  assess  these  nonisothermal  effects.  However,  we 
remind  the  reader  that,  under  thermophysical  conditions  such 
that  Bit,«  1,  then  our  isothermal  -results  (Figures  7  and  8) 
could  be  taken  over  to  find  even  under  nonisothermal 
conditions,  since  the  dominant  temperature  nonuniformities 
would  then  occur  only  in  the  external  (thermal)  gaseous  bound¬ 
ary  layer. 

Conclusions 

Principal  conclusions  for  Dj<3  aggregates 

A  simple  pseudo-continuum  method  has  been  developed, 
described  and  illustrated  for  calculating  the  accessible  surface 
area  of  a  large  open  aggregate  containing  N  primary  spherical 
particles  (A^»l)  suspended  in  a  gas  at  pressures  such  that 
the  mean-free-path  is  not  necessarily  negligible  or  appreciable 
compared  to  the  overall  aggregate  radius.  For  this  purpose, 
the  A/-pariicIe  aggregate  has  been  treated  as  a  porous  granular 
solid  characterized  by  the  scaling  exponent  Dj-  dlnN/dlnr, 
and  the  primary  sphere  diameter  2/?,.  A  simple  yet  rational 
correlation  scheme  has  been  developed,  based  on  the  more 
familiar  Z7y=3  case  and  the  notion  of  an  effective  (stretched) 
Thiele  modulus  which  can  predict  to  within  about  3%  the 
numerically  calculated  normalized  accessible  surface  area  of 
such  an  aggregate  over  the  entire  range  of  parameters  of  phys¬ 
ical  interest,  viz.:  primary  panicle  number,  N\  structural  ex¬ 
ponent,  Df,  reaction  probability,  a;  and  the  ambient 

gas  mean-free-path  (measured  in  primary  particle  diameters). 

As  a  byproduct,  valuable  insight  has  been  provided  on  recent 
experimental  techniques  for  in  situ  probing  of  the  accessible 
area  of  fractal  aggregates,  for  example,  the  scavenging  of  ra¬ 
dioactive  lead  atoms  and  photoelectron  emission. 


Application  of  correlations  to  a"l,  calculations  for  co¬ 
agulation-aged  suspensions 

This  formulation  and  the  above-mentioned  correlation  open 
the  door  to  calculations  of  the  accessible  area  per  unit  volume 
for  coagulation-aged  populations  of  aggregates,  as  observed 
in  many  technologies,  and  in  natural  environments.  Thus,  if 
there  are  a  total  of  aggregates  per  unit  volume,  with  spec¬ 
ified  normalized  size  distribution  pdfiN} ,  we  can  now  rapidly 
calculate  (at  least  for,  say,  /V>30)  the  most  difficult  part  of 
the  integrand  appearing  in: 

=  I  fli'V  >)ovc,.M-(A';  Dy,  Of,  pdfiN^dN  (49) 

and  compare  the  result  of  this  integration  to  the  total  area  per 
unit  volume,  associated  with  the  sum  of  all  primary  par¬ 
ticles  present  irrespective  of  their  state  of  aggregation.  Indeed, 
this  -ratio,  which  can  be  considered  an  overall  effec¬ 

tiveness  factor  for  the  aggregate  population,  will  be  seen  to 
dictate  aggregated  particle/carrier  fluid  exchange  processes  in 
many  applications  involving  flocculated  suspensions  (Rosner 
and  Tandon,  1993a),  including  physical-  and  chemical-vapor 
deposition  (Castillo  and  Rosner,  1988;  Rosner  and  Liang,  1987) 
and  “particle-modified”  chemical  vapor  deposition  (Komi- 
yama  et  al.,  1991). 

Future  work  on  transport  to/from  large  aggregates 

A  number  of  the  generalizations  outlined  in  the  fourth 
section  may  prove  to  be  necessary  and  ultimately  valuable, 
building  on  the  effective  porous  object  results  and  approach 
outlined  here.  Indeed,  to  belter  understand  the  large  N-limit 
this  approach  is  being  extended  to  all  interesting  mass-,  mo¬ 
mentum-,  and  energy  exchange  properties  of  aggregates,  in¬ 
cluding  their  ihermophoretic  “propulsion”  (Rosner  et  al.,  1991) 
and  radiative  properties  (Dobbins  and  Megaridis,  1991;  Mac¬ 
kowski,  1992b).  It  would  also  be  interesting  to  investigate  the 
vapor  nucleating  ability  of  such  aggregates,  including  uptake 
rates  by  capillary  condensation  in  the  singular  “pores”  created 
where  primary  particles  touch.  These  investigations  of  aggre¬ 
gates  containing  many  primary  spherical  particles,  including 
their  rearrangement  (restructuring)  kinetics  (Cohen  and  Ros¬ 
ner,  1994)  are  now  underway,  using  some  of  the  results  of  this 
article  and  analogous  pseudo-continuum  methods  (Rosner, 
Cohen,  and  Tandon,  1993). 
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Notation 

Oi  -  surface  area  of  one  primary  sphere  (Figure  2) 
o'”  =  surface  area  per  unit  volume 
A  =  reactant  species  lOH(g),  C2H:(g),  ^"Pb(g)),  or  area 
M  =  aspect  ratio  (of  prolate  spheroid) 
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Bi, 


c 

Ca 

D 


D, 

f 


^eff 

k” 

Kn 

I 

s 


/?, 

Bt,^r 

^max 

V 

X 


Biot  number  for  heat-transfer  (ratio  of  internal*  to  ex¬ 
ternal  heat-transfer  resistance) 
normalized  species  A  local  concentration  variable 
mean  thermal  speed  of  molecules  A,  [^kgT/ {rm 
normalized  effective  diffusion  coefficient;  Eqs.  2~18 
effective  diffusion  coefficient  for  species  A  in  local  po¬ 
rous  medium 

coefficient  for  diffusion  of  species  A  through  prevailing 
vapor  mixture 

fractal  dimension,  exponent  dlnN/dlnr 

Knudsen  diffusion  coefficient  (single  straight  pore  of 

average  diameter) 

fraction  of  incident  molecules  scattered  diffusely  (by  sur¬ 
face) 

effective  pseudo-homogeneous  first-order  rate  constant 

first-order  heterogeneous  rate  constant 
Boltzmann  constant 

Knudsen  number  (ratio  of  gas  mean-free-path  to  char¬ 
acteristic  dimension) 
gas  mean-free  path 
molecular  mass  of  species  A 
local  number  of  density  of  species  A 
total  number  of  primary  spheres  in  aggregate  (see  Figure 
2) 

Nussell  (Sherwood)  number  for  mass  transfer  to/from 
sphere 

radius  (measured  from  aggregate  center-of-mass);  Figure 
2 

radius  of  primary  spheres;  Figure  2 
inner  core  radius;  Eq.  3 
radius  of  gyration;  Figure  2 
maximum  (outer)  radius  defined  by  Eq.  4 
absolute  temperature  of  solid  at  r= 
volume 

argument  (of  the  indicated  function): 


Greek  letters 


a 


0 


r 

V 


X 


reaction  probability,  condensation  coefficient,  and  so  on 
(0<a<l) 

filling  factor  (Eq  =  1  -  tum 

local  void  fraction  (1  -d>) 
normalized  radial  position  variable, 

“effectiveness”  factor:  accessible  area  as  a  fraction  of 
Na, 

function  defined  by  Eq.  31 
radiation  wavelength 

tortuosity  of  local  pore  structure  defined  by  Eqs.  6  and 
7 

Thiele  modulus  (reciprocal  of  dimensionless  reactant 
penetration  depth) 
solid  fraction  (1  -  e) 


Miscellaneous 

BC  =  boundary  condition 
cm  =  center  of  mass;  Figure  2 
FAM  =  finite-analytic  method 
ODE  =  ordinary  differential  equation 
0-(  >  =  order  of  magnitude  “operator” 

S-0  =  Schmidt-Ott 
tanh-(  >  =  hyperbolic  tangent  function 

=  radioactive  isotope  of  the  element  lead 
(~)  =  normalized  by  conditions  at  r= /?  mat 
<  >5  =  average  value  defined  by  Eq.  33 
O  =  argument  of  indicated  function 
fet-i  )■  =  function  of  indicated  argument 

pdfi  >  =  probability  density  function  of  indicated  random  vari¬ 
able 


Subscripts 

OSS  =  pertaining  to  aggregate 


A 

c 

eff 

eq 

exp 

ext 

f 

Syr 

h 

int 

K 

lim 

m 

mix 

max 

pore 

Q 

rot 

2/?, 

I 


pertaining  to  species  A 

pertaining  to  core  or  to  continuum  (high-pressure  limit) 

conditions 

effective 

pertaining  to  local  equilibrium 

experimental 

external 

fractal 

gyration 

pertaining  to  heat  transfer 

internal 

Knudsen 

limiting  value 

pertaining  to  mass  transfer 

pertaining  to  local  vapor  mixture 

maximum  value 

based  on  (mean)  pore  diameter 

index  appearing  in  Eq.  33;  ^^>0 

rotational 

based  on  diameter  of  primary  particles 

primary  particle  (spherule);  Figure  2 

evaluated  at  r=oo  (far  from  aggregate  center-of-mass) 
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SUMMARY  - 

We  review  here  results  of  our  recent  research  on  the  deposition  dynamics  of  combustion- generated  particles  in 
power  production  and  materials  synthesis/processing  technologies.  In  this  brief  summary,  which  emphasizes  the 
capture  of  inertially  impacting  particles*,  we  outline  and  illustrate  the  results  of  recently  developed  methods  to  predict 
total  surface  deposition  rates  and  associated  convective  heat  transfer  reductions  for  targets  exposed  to  a  distribution  of 
particles  suspended  in  a  mainstream.  Our  methods  combine  the  essential  features  of  recently  developed  single  particle 
sticking  probability  laws  with  correlations  of  the  inertial  impaction  of  particles  on  targets  in  high  Reynolds  number 
cross-flow,  to  develop  formulae  and  ’universal'  graphs  which  provide  the  dependence  of  particle  deposition  rates,  and 
associated  reductions  in  convective  heat  transfer,  on  such  system  parameters  as  mainstream  velocity,  mean  suspended 
particle  size  and  target  size.  The  deposition  rate/deposit  microstructue/properties  prediction  and  correlation  procedures 
illustrated  here  [R  9,10]  can  be  used  in  many  ways;  eg,  incorporated  into  improved  'fouling  propensity  indices’,  to 
motivate,  evaluate  and  implement  "ruggedization”  and/or  fouling  reduction  strategies,  and/or  to  optimize  the 
"harvesting”  of  tailored  particles  for  subsequent  processing. 

Details  on  our  theoretical  studies,  and  their  immediate  laboratory  antecedents  [L1,R4,5,16],  will  be  found  in 
our  archival  references  (Section  8).  A  judicious  blend  of  (numerical  and  physical)  experiments,  theory,  and  intuition  will 
continue  to  be  necessary  to  economically  arrive  at  methods/results  to  improve  the  generality  and  accuracy  of  particle- 
deposition-related  design  calculations  for  a  wide  variety  of  equipment/fuel  types.  Our  computational  and  correlation 
methods  are  being  extended  to  treat  more  complex  situations  of  practical  importance  in  power  generation  and  materials 
synthesis/processing  applications. 


1.  INTRODUCTION.  BACKGROUND 

Initially  clean  heat  exchanger,  containment,  or  target  surfaces  exposed  to  high  temperature 
flowing  suspensions— eg.,  ash  or  soot  particles  in  fossil  fuel  (oil  or  coal-)  combustion  products,  can 
acquire  a  sufficient  fraction  of  this  solid  material  to  cause  a  noticeable  decline  in  heat  transfer  or 
aerodynamic  performance.  This  decline  is  associated  with  the  local  growth  of  microgranular 
insulating  layers,  often  undesired  and  creating  the  need  for  periodic  shutdown  for  cleaning.  To 
answer  such  questions  as:  required  maintainence  intervals  for  a  particular  (class  of)  fuel(s) ;  assess  the 
most  economical  degree  of  fuel  'cleaning';  or,  in  particle  synthesis/processing  applications,  to  select 
optimum  conditions  for  harvesting  desirable  particles  or  coatings,...  a  quantitative  understanding  of 
deposition  rates  /deposit  microstructure/properties  is  clearly  necessary.  From  the  viewpoint  of 
capturing  a  non-negligible  fraction  of  the  mainstream  particle  flow  rate,  inertial  impaction  is  usually 
the  mechanism  responsible  for  most  of  the  deposited  mass  and  volume.  However,  by  far  the  greatest 
uncertainty  in  making  such  deposition  rate  predictions  is  associated  with  the  appropriate  single 

impacting  particle  capture  fraction  ,  or  sticking  probability  S  [RIO].  In  general,  5  is  a  function  of 
both  incident  velocity  and  angle,  but  not  yet  fully  understood  even  for  particle  impaction  on  'dry' 

granular  deposits  [K3,  R  10  ].  Moreover,  in  most  applications  s  will  be  determined  in  part  by  the 
simultaneous  rates  of  vapor  and  much  smaller  particle  (non-inertial)  deposition  [R6]*.  As  summarized 
below,  we  have  extended/applied  inertial  impaction  theory  [FI,  I1,K2,  R13]  and  applied  what  we 
now  know  about  single  particle  capture  behavior  [  K3  ,  RIO  ]  to  predict  the  sensitivity  of  deposition 
rate  to  system  parameters  {eg.,  mainstream  velocity,  particle  mass  loading,  rebound  velocity,  mean 


*For  a  review  of  multi-component  vapox  convective-diffusion  and  thermophoretic  particle  mass  transfer  in  chemically  reactive 
flow'  systems,  see,  e.g.,  [R8,15,17  ].  Our  prior  experimental  and  theoretical  studies  of  alkali  vapor  deposition  are  conveniently 
summarized  in  [R  4  ].  For  a  useful  overview  of  heat  exchanger  fouling  research  as  of  ca.  1989  see:  [M  1  ]. 


particle  size,  target  diameter,  etc.).  This  capability  can  clearly  be  used  to  motivate,  evaluate  and 
implement  "ruggedization"  and/or  fouling  reduction  strategies,  or  optimize  the  'harvesting'  of  tailored 
particles  for  subsequent  processing. 

2.  PREDICTION  of  PARTICLE  DEPOSITION  RATES 

We  have  developed  a  convenient  formalism  [R9]  for  making  rational  engineering  predictions 
of  deposition  rates  in  high-gas  velocity  particle-laden  environments,  based,  in  part,  on  recently 
developed  single  particle  sticking  capture  laws  [K3,R10]  for  impaction  on  granular  deposits. 
Exploiting  such  information,  even  when  available  for  the  particular  materials  combinations  of  interest, 
to  anticipate  deposition  rates  in  realistic  engineering  environments  was  previously  thought  to  be  a 
computationally  demanding  task  since  it  is  necessary  to  track  the  impingement  frequency,  velocity  and 
incidence  angles  of  all  the  different  size  particles  in  the  mainstream  capable  of  striking  the  target 
locations  of  interest,  invoking  the  abovementioned  sticking  ( or  'rebound')  laws  at  each  such  point  to 
predict  the  corresponding  cumulative  local  deposition  rate.  However,  by  focusing  our  attention  on 
specific  canonical  geometries  {eg.,  circular-cylindrical  targets  in  cross-flow  (Fig.l))  and  introducing 
a  modest  number  of  defensible  approximations  (Section  3)  to  summarize  the  predicted  and/or 
measured  rebound  behavior  on  the  solid  surfaces  of  interest,  we  have  shown  that  the  tedious  portion 
of  such  deposition  rate  predictions  can  be  efficiently  carried  out  'once-and-for-all',  thereby  reducing 
the  engineering  problem  of  predicting  target  deposition  rates,  and  the  associated  reduction  in 
convective  heat  transfer  rate,  to  that  of  multiplying  a  readily  calculated  reference  deposition  rate  by  a 
set  of  'universal'  dimensionless  deposition  rate  functions,  some  of  which  are  illustrated  here.  Our 
reference  deposition  rate  is  that  which  would  be  expected  in  the  prevailing  environment  if  all 
mainstream  particles  had  the  mean  size  and  were  captured  upon  impacting  the  target  surface  without 
experiencing  aerodynamic  deceleration  or  deflection  due  to  the  carrier  gas  flow.  For  convenience  and 
generality,  our  results  are  cast  in  terms  of  the  following  dimensionless  parameters:  ratio  of 
mainstream  velocity  U  to  the  reference  critical  velocity,  Vp_crit(  '^crit  )>  for  particle  rebound  from  the 
solid  surface,  ratio  of  mean  abrasive  particle  size,  v,  to  the  threshold  ('critical')  size,  Vcrit^  required 

for  impaction  on  the  target  in  the  prevailing  flow  environment  [R3],  spread,  CTg  ,  of  the  mainstream 
particle  size  distribution  Coofv)  (here  assumed  'log-normar[R3]),  and  the  characteristic  'slip 
Reynolds  number  for  the  critical  size  particles  in  the  mainstream  [II].  In  this  way  we  have  shown 
that  many  previously  observed  characteristics  of  granular  deposits  (eg.,  fouling  layers  layers) 
including  their  frequent  "lobular"  appearance  (Fig.  3),  can  be  understood  theoretically. 

3.  PRINCIPAL  ASSUMPTIONS 

To  incorporate  the  essential  physical  phenomena  without  making  unrealistic  idealizations, 
most  work  to  date  has  been  based  on  the  following  basic  assumptions: 

At  Local  particle  impaction  frequencies,  velocities,  and  angle-of-incidence  can  be  calculated  with  sufficient 
accuracy  from  recent  correlations  summarizing  the  results  of  individual  su.spended  non-Brownian  particle  trajectories 
calculated  for  steady,  in  viscid  flow  past  the  target ,  including  non-Stokesian  drag  corrections  (c/  [I  l](Fig.  1) 

A2  Even  for  impaction  on  granular  deposits  [RIO],  single  particle  capture  probability  laws  at  particular 
velocities  and  incidence  angles  can  be  invoked  to  predict  average  deposition  rates  in  engineering  applications  where 
suspended  particles  of  different  sizes  arrive  over  a  broad  range  of  impact  velocities  and  incidence  angles. 

A3.  ’Rebounding’  particles  donot  appreciably  influence  incoming  particles,  nor  deposit  in  appreciable  numbers  upon 
re-impaction  on  the  same  target 

A4  Predicted  'initial’  deposition  rate  trends  (spatial  distributions)  determined  on  an  initially  smooth  target  can  be 
used  to  anticipate  longer  time  deposition  rate  trends  on  inevitably  roughened  targets  which  ultimately  depart  from  their 
original  shape  due  to  localized  deposit  growth. 

A5  The  mainstream  population  of  suspended  particles  is  approximately  log-normal  with  respect  to  particle  volume 
[R3]  and,  while  the  particle  mass  loading,  cOp,  in  the  mainstreamlneed  not  be  very  small  [R1  ,  PI  ],  the  volume 
fraction,  (|)p,  corresponding  to  the  total  particle  number  density  Np  and  mean  particle  volume  v  is  negligible. 


Subject  to  these  assumptions,  as  noted  above  actual  local  deposition  rates  at  any  position  x/Lt 
(=polar  angle  6  for  a  circular  cylinder  (Fig  1)  if  Lt  is  taken  to  be  dt/2)  can  be  expressed  as  the 
product  of  an  easily  calculated  reference  deposition  rate,  (DR)|.gf  and  a  universal  dimensionless 
function  D(x/Lt,...)  introduced,  calculated(Section  4 )  and  plotted  (Fig.  3)  over  the  interesting  range 
of  mean  suspended  particle  diameters  (expressed  as  a  multiple  of  critic^  diameter  required  for  inertial 
impaction  in  the  prevailing  environment  [R3])  and  a  dimensionless  velocity  ratio  characterizing  the 
rebound  behavior  of  the  particles  on  the  surface  of  interest.  An  attractive  feature  of  this  type  of 
formulation  is  that  the  dimensionless  function  D,  and  its  average  value  D  over  the  'upwind-facing' 
surface  of  the  target,  can  be  calculated  'once-and-for-all'  via  straightforward  numerical  quadratures  in 
terms  of  an  acceptably  small  number  of  dimensionless  parameters  defining  the  system  (application). 
Indeed,  the  availability  of  such  results  dramatically  simplifies  the  task  of  predicting  local  and  total 
deposition  rates  for,  say,  heat  exchanger  tubes  in  the  cross-flow  of  ash-laden  combustion  products 
[R9] .  Moreover,  by  introducing  a  modified  (heat  transfer  coefficient-weighted-)  average  value  of  D, 
designated  (Eq.2(below),Fig.4),the  prediction  of  the  associated  convective  heat  transfer  reduction. 

has  also  been  dramatically  simplified  [R9]. 


4  SINGLE  PARTICLE  CAPTURE  PROBABILITY  LAWS 
and  the  CALCULATION  of  D(x/Lt,  ...) 

We  have  applied  the  micromechanical  theo^  of  particle  capture  sticking  fraction  (which 
provides  the  functional  form  of  s  when  particular  projectile  particles  are  directed  at  particular  target 

materials  (including  granular  deposits)  at  a  known  velocity  Vp  and  angle  of  incidence  6;  (cf  the  target 
outward  normal)  to  predict  local  and  total  particle  capture  rates  for  specific  cases  of  practical 
importance;  eg.  a  cylindrical  target  immersed  in  a  polydispersed  suspension  of  such  particles [R3, 9  ]. 
Three  distinct  classes  of  single  particle  capture  laws  have  been  considered,  as  follows:  constant 
capture  fraction,  "on-off  (Fig. 2a)  capture  behavior  (expected  on  a  'clean'  (particle-free)  smoothe 
surface),  and  capture  by  a  'granular'  deposit  (above  the  particle  critical  velocity,  the  sticking 
probability  does  not  fail  abruptly  to  zero  but,  rather,  exhibits  an  exponential  'tail'  ([K3,R10]  see  Fig. 
2b). 

The  dimensionless  local  deposition  rate  function  D(x/Lt)  is  explicitly  given  by  an  integral  of 
the  form: 

D(x/Lt,...)  =  i*Js(v,x/L,)Tliocai(v,x/Lt)v-C«(v)-dv  (O 

0 

To  complete  the  calculation  of  D(x/Lt,...>in  any  particular  situation  (see,  eg.,  Fig.3)  we  inust  specify 
the  three  inertial  impaction  functions:  impact  velocity ,Vp(v,  x/Lt)/U>  angle-of-incidence  (Fig.  1) 

0j<v,  x/Lt>  (these  determine  the  local  capture  fraction)  and  the  local  dimensionless  particle 

impingement  frequency  primary  goal  has  often  been  the  mean  value  D  of  D  over 

the  upwind-facing  surface  of  the  target.  For  a  circular  cylinder  it  is  interesting  to  note  that  D  so 
defined  is  mere!;  2/n  times  the  conventionally  defined  total  capture  (efficiency)  fraction,  based 

on  target  'frontal  area'  [R15]. 


5.  ASSOCIATED  REDTICTION  in  CONVECTIVE  HEAT  TRANSFER  RATE 

One  frequent  reason  for  interest  in  deposition  dynamics  is  the  consequence  of  even  a  thin 
deposit  for  heat  transfer  performance.  Accordingly,  we  derived  the  following  interesting  explicit 
relation  between  fractional  convective  heat  transfer  reduction  and  the  abovementioned  local 
deposition  rate  function.  An  important  result  of  our  analysis  is  the  relevance  of  the  (Nuj^fx/Lt)- 
'weighted'-)  integral  : 
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where  the  integrand  contains  the  position  dependence  of  the  local  convective  heat  transfer  coefficient 
and  xm  is  the  maximum  value  of  x  capable  of  experiencing  particle  impacts  for  a  symmetrical  target 
The  relevance  of  Du  to  heat  transfer  reduction  is  clearly  understood  when  it  is  recognized  that  equal 
thickness  local  deposits  will  have  greater  effect  on  overall  convective  heat  transfer  if  located  near  the 
position  of  maximum  Nu^;  eg.  the  forward  stagnation  point  region  (where  NUh<e;Re)  maximizes  in 
the  Re-range  of  interest)  than  if  located  near  0=71/2  radians  (c/  Figs.  1,3).  The  fractional  reduction  of 
convective  heat  transfer  rate  with  time  t  on  stream  is  found  to  be  given  by  [R9]. 


l(DR)ref 


('’»), .0 


dt 


-Aq'w  -  1 


^dep 


-NUh{0) 


Nu. 


(3) 


To  summarize,  for  the  canonical  case  of  say,  a  circular  cylinder  in  Retl/2»1  crossflow,  our 
'universal'  results  for  D(0>and  ,  coupled  with  Eq.(3  )  and  a  rational  estimate  of  (l-<e>).kjep 
appropriate  to  the  granular  deposit  in  question,  can  be  used  to  predict  the  time-on-stream 
con-esponding  to,  say,  a  specified  fractional  loss  of  convective  heat  transfer  rate  of  9w 

(eg.,10%).  Under  conditions  for  which  the  deposits  are  not  necessarily  confined  to  the  forward 
stagnation  region  {eg.,  for  impact/rebound  on  a  bar^surface  at  large  values  of  (v/v^^lj)  an 
U/V  ),  then  Dh  is  found  to  be  noticeably  smaller  than  D.  .  •  j 

^’^"^Figure  3  is  a  representative  polar  plot  of  dimensionless  deposition  rates  on  the  opwin 
surfaces  of  the  cylinder  in  crossflow  for  one  of  the  sticking  laws  considered  (Fig.  2).  The  mbound 
parameter  (U/Vj,^(Vj,^»  is  only  0. 1  for  the  case  (shown)  of  capture  on  clean  solid  surfaces  . 

calculated  contours  are  for  different  values  of  dimensionless  mean  volunie,  |,  defined  2S  the  ratio  of 
V  of  the  mainstream  distribution  and  the  critical  particle  volume,  v^^it-  For  the  case  of  clean  s 
surfaces,  an  impacting  particle  is  captured  only  if  the  normal  velocity  component,  Vp  is  ess  m 
the  size-dependent  critical  velocity,  V  „i,(v>.  Therefore,  as  the  dimensionless  mean  size  of  the 
particles  increases,  a  larger  fraction  of  these  particles  rebound  in  the  forward  stagnation  region  but 
somewhat  smaller  particles  are  captured  at  larger  angles,  0,  leading  to  the  formation  of  deposits  with 

nff-axis 'lobes' (Fie.  3),  which /lave  been  frequently  observed.  _ 

FigL  4  suminarizes  representative  results  for  the  function  Dh  needed  ^ 

heat  transfer  reductions  {via  Eq(  3  )).  trends  for  Dh  are  the  same  as  those  for  D  but,  as  expected, 
there  is  a  systematic  decrease  below  D  in  those  cases  for  which  there  is  appreciable  off-axis 

byproduct  of  this  approach  is  the  ability  to  evaluate  the  local  sensitivity  of 
deposition  rates  to  key  system  parameters  [R9],  such  as  gas  velocity,  mean  particle  size  in  the 
mainstream  and  target  dimension.  Our  results  indicate  the  presence  of  interesting  opposing 

tendencies _ for  example,  an  increase  in  gas  velocity  increases  the  particle  impingement  rate  but  also 

increases  the  fraction  of  impacts  leading  to  rebound.  Sirnilarly,  an  increase  in  target  dimension  wou 
reduce  the  frequency  of  impacts  but  increase  the  likelihood  of  capture  upon  each  (lower  y^e  oci  y) 
[Sl^act.  ae^y'^  of  these  'tradeoffs'  is  necessary  to  evaluate,  say,  the  efficacy  of 

rival  strategies  for  reducing  heat  exchanger  fouling  rates. 

6  nFPO.SlT  MTCROSTRIJCTURE/THERMOPHYSICAL  PROPERTIES 

*  In  many  applications  it  is  necessary  to  embrace  phenomena  explicitly  trussing  from  the 
■baseline'  analyL  outlined  above  {eg.,  larger  particle 

nresence  and  capture  of  condensible  vapor  and/or  submicron  (subcntica  )  particles  fsee,  eg., 
rri  R6  101).  For  example  we  would  expect  systematic  shifts  in  the  critical  velocity  for  particle 
rebound,  V  relative  amount  of  vapor  and/or  sub-micron  (supercntical)  particle 

deposition,  locally. 
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Additionally,  it  is  clearly  necessary  to  understand  the  interesting  links  bet\  on  deposit 
microstructure/properties  and  particle  deposition  mechanism  [Tl,5].  Thus,  in  comparing  such 
deposition  rate  predictions  with  experiments,  it  should  be  borne  in  mind  that,  for  'thick'  granular 
deposits,  one  must  incorporate  information  on  the  microstructure  including,  at  least,  the  void 

fraction  <e>of  such  deposits  (see,  eg.,[Tl,4,R10]).  We  know  from  these  ancillary  studies  that  the 
local  deposit  thickness  and  associated  local  thermal  resistance  at  any  given  time  on-stream  are  not 
determined  only  by  the  history  of  the  local  particle  volume  fluxes  due  to  impaction/capture. 

The  particle-level  dynamical  simulations  which  led  to  our  currently  used  capture  fraction  laws 
for  particle  impact  on  dry  mono-sized  deposits  (Section  4,  Fig.  2,  and  [Kl,  RIO])  also  lead  to 
particular  microstructures  [RIO]  which  have  been  characterized  using  a  number  of  distribution 
functions  [T4  ]  including  'pore'-size  and  contact  normals.  These  statistic^  descriptors  have  also  been 
used  for  deposits  grown  using  simpler  'algorithms'  regarding  arrival  mechanism  and  'rolling  events' 
prior  to  incident  particle  'arrest'  (Fig.  5,  after  [Tl-5]).  Using  accelerated  random-walk  (Brownian 
mass-  or  heat-  'tracer')  simulation  methods  we  have  also  developed  and  illustrated  methods  to  predict 
the  anisotropic  effective  'diffusivities'  of  these  deposits,  including  the  effective  thermal  conductivity 
kdep,eff  (Fig-  6,  after  [T3])  and  the  effective  Fick  diffusivity  for  vapor  diffusion  through  the  gas- 
filled  interstitial  spaces  between  the  spheres.  We  are  currently  generalizing  these  methods  to  include 
the  effects  of  particle  'polydispersity',  capillary  condensation,  and  partial  sintering  ('consolidation'). 

7.  CONCLUSIONS 

Despite  the  daunting  complexities  remaining  to  be  overcome  in  the  design  and  operation  of  power  plants  (eg., 
utilizing  a  broad  spectrum  of  ash-bearing  fuels,  including  refuse-based  ’fuels'),  as  well  as  combustion-derived  particle 
synthesis/processing  devices  {eg.,  optical  wave  guide  manufacture  [Rl]),  these  recent  methods  and  results  are  indicative 
of  the  potentially  useful  simplifications  and  generalizations  emerging  from  current  fundamental  research  studies  of 
suspended  particle  deposition  mechanisms  and  their  connection  to  microparticulate  deposit  formation.  Long-range 
investigations  should  continue  in  each  of  the  underlying  theoretical  areas  summarized  above,  as  well  as  others  beyond 
the  scope  of  this  brief  paper  (see,eg.,[C4,Gl,2,  R7,8,l  1,13,14],  Only  then  will  it  be  possible  to  provide  engineers 
with  significantly  better  "tools"  for  making  rational  assessments  of  particulate  deposition  phenomena  in  future,  high- 
performance  equipment.  We  hope  this  brief  overview  and  the  reasonably  complete  bibliography  of  our  work 
(including  titles;  Section  8)  will  accelerate  this  process  via  the  requisite  and  timely  exchange  of  R&D  information 
among  PTF-participants. 
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Fig.  1  Fouling  of  a  heat  exchanger  surface  (here  circular 
cylinder  in  high  Reynolds  number  crossflow)  in  response 
to  the  arrival  of  impacting  particles  log-normally 
distributed  with  respect  to  size.  Particle  capture  fraction  s 
dependent  on  incident  velocity  (Fig.  2)  and  angle  Gj  [R  9  J 


a  circular  cylinder  in  crossflow  for  several  values  of  the 
particle  population  size  parameter  v/v’eriti  Ca.ses  shown: 
Threshold  velocity  for  critical  size  particle  rebound  equal 
to  ten  times  the  mainstream  velocity,  U  (after  [R91).Nole 
tendency  to  form  off-axis  'lobes’  when  v/verit  exceeds 
about  100. 


Fig. 5  Particle-level  computer  simulation  (2D)  of 
microstructures  for  granular  deposits  of  uniform  size 
spoherical  particles  (after  [TI.4]).  Cases  ahown: 
algorithmic  deposits  with  various  degrees  of 
restructuring"  (0-,  1 -.2-,3-"roIling"  events  after  initial 
impact;  normal  ballistic' incidence) 


Fig.  2  Single  particle  capture  probability  laws  (particle 
velocity  dependence)  considered  |R9|.  Cases  shown: 
capture  on  a  clean  solid  surface;  and  incident  particle 
capture  by  a  granular  deposit  (after  IK31) 


Fig.  4  Upwind  cylinder  ,<;urf,nce-aver.nged  deposition  mle 
lunclion  dictating  convective  heat  transfer  reductions  {via 

Eq.  (3  )).  Dependence  on  mean  particle  size  (volume)  in 
the  mainstream  population  for  each  class  of  particle 
sticktng  fraction  law  (Fig.  2).  (after  [R  9)) 


Fig,  6  Effective  orientation-averaged  thermal  conductivity 
of  anisotropic  computer-simulated  granular  deposit's 
comprised  of  uniform-sized  .spherical  particles  (alter  Refs. 
[Tl,3)).  Brownian  'heat  tracer'  simulation  results  as  a 
function  of  'void '/particle  thermal  conductivity  ratio.  Also 
shown:  experimental  data  and  rigorous  theoretical  bounds. 
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